TIDAL  FLOW  IN  ENTRANCES;  WATER-LEVEL 
FLUCTUATIONS  OF  BASINS  IN  COMMUNICATION 
WITH  SEAS 

G .  H .  Keulegan 

Committee  on  Tidal  Hydraulics  (Army) 
Washington,  D.  C. 

July  1967 


DISTRIBUTED  BY: 


Nattenal  Technical  Infoneatjoc  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


riVHNICAl  OUi  u  riN  NO.  14 


TIDAL  FLOW  IN  ENTRANCES 

WATER-LEVEL  FLUCTUATIONS  OF  BASINS 
IN  COMMUNICATION  WITH  SEAS 

G.  H.  Keultgan 


July  1967 


NATIONAL  TECHNICAL 
INI  oi -MAT ION  SERVICE 

■.  i  .  -  • ' '  ■  ' 


D  D  C. 

Li; 


u 


\m 


lliiilatanLi’lS! 

- . a  ^ 


Committo®  on  Tidal  Hydraulics 
CORPS  OF  ENGINEERS,  U.  S.  ARMY 


DISTRIBUTION  OF  THIS  DOCUMENT  IS  UNLIMITED 


Best  Available  Copy 


1 1  0H  hdal  MroRAUUcs 


|M>#I  N  -■>. 

TltU 

Dote 

Fvaluulitn  of  Pressnt  State  of  Km  pledge  of  t  art«*  Affec'ing  Tidal  Hydroullet 
and  Related  Phenomena 

r.b.  1950 

'i 

Bibliography  on  Tidal  Hydraulic* 

*eb.  1954 

Supplement  Ho.  1,  Material  Compiled  Through  May  1955 

June  1955 

Supplement  No.  7,  Material  Compiled  from  May  1955  to  May  1957 

fAoy  1957 

Supplement  No.  3,  Material  Csmplled  frost  May  1957  to  May  1959 

Moy  ’,959 

Supplement  No.  4,  Material  Compiled  from  May  1959  to  May  1965 

May  1965 

? 

Evaluation  of  Present  State  of  Knowledge  of  Ffletors  Affecting  Tidal  Hydraulics 
a<td  Related  Phenomena  (revised  edition  ol  Report  No,  1) 

Moy  1965 

T  •chiiicc.  I 

Bulletin  No.  1  Itld 


I 

7 

3 

4 

5 

6 

7 

8 
9 

10 

11 


ir 

13 

14 


"Vi*.  ¥4>-'  -I- 


M 


Sed  ment  Discharge  Measurements  in  T idol  Waterway* 

Fr>i»h  Water-Salt  Water  Density  Currontf,,  o  M o\ur  Causa  of  Slltotlon  In  Estuaries 
Tidal  Flow  In  Entrances 

Soil  a*  a  Factor  in  Shoaling  Proe**«.<is,  a  Literature  Revinw 

Ore-Olmentlonal  Analysis  of  SrJinlfy  Intrusion  In  Estuor'ins 

Typical  Major  Tidal  Hydraulic  Problem*  In  United  Staten  and  Research  Sponsored 
by  the  Corps  of  Engineer*  Committee  on  Tidal  Hydraulic* 

/.  Study  of  Rheologi*  Properties  ol  Estuarial  Sediments 

Chonnel  Depth  cs  a  Factor  In  Estuarine  Sedimentation 

A  Comparison  of  on  Ettuory  Tide  Calculation  by  Hydioullr,  Model  an  I  Computor 
Significance  nf  Clay  Mineral*  In  Shoaling  Problem* 

Extract*  from  The  Manual  of  Tide* 

’Unpublished  Consultation  Report*  on  Carp*  of  Engineer*  Tidal  Pro|ect* 

Two-Dimensional  Arpect*  of  Salinity  Intrusion  in  Estuorles:  Analysis  of  Salinity 
ond  Velocity  Distribution* 

Tldol  Flow  In  Entrances;  Woter-Lovel  Fluctuation*  ol  Bosln*  In  Communication 
with  Seas 


°^1 


n 
1 1 


■J 


■ft 


Destroy  this  raport  wh«rt  no  longer  needed.  Do  not  return 
it  to  the  originator. 


Dotr 

May  1954 
April  1957 
Jon.  1960 
Jun#  I960 
June  1965 

June  1963 
Sept.  1963 
Mar.  1965 
June  1965 
Sept.  1966 
Sept.  1666 
Dec.  I 

June  19’ 7 

July  V-  f 


I  The  finding*  *n  this  report  ore  not  to  be  construed  o*  on  officio! 
Deportment  of  the  Army  position  unless  »o  deelgnoted 
by  other  outhorlred  document*. 


Best  Available  Copy 


TECHNICAL  BULLETIN  NO.  14 


TIDAL  FLOW  IN  ENTRANCES 

WATER-LEVEL  FLUCTUATIONS  OF  BASINS 
IN  COMMUNICATION  WITH  SEAS 

by 

G.  H.  Keulegan 


Committee  on  Tidal  Hydraulics 
CORPS  OF  ENGINEERS,  U.  S.  ARMY 

IC  VICKtIUftO.  Mlt« 


DISTRIBUTION  OF  THIS  DOCUMENT  IS  UNLIMITED 

j 


PRESENT  MEMBERSHIP  OF 
COMMITTEE  ON  TIDAL  HYDRAULICS 


Members 

J.  B.  Tiffany 
J.  B.  Lockett 
P.  A.  Becnel,  Jr, 

J.  M.  Caldwell 
E.  B.  Conner 
A.  B.  Davis,  Jr. 

J.  W,  Harris 
J.  B.  McAleer 
J.  F.  Phillips 
T.  Saville,  Jr, 

E.  A.  Schultz 
H.  B.  Simmons 

Liaison 
J.  H.  Douma 

Consultants 

Dr.  Arthur  T.  Ippen 
Dr.  G.  H.  Keulegan 
Dr.  D.  W.  Pritchard 
Mr.  C.  F.  Wicker 


Waterways  Experiment  Station  Chairman 

North  Pacific  Division  Recorder 

New  Orleans  District 

Coastal  Engineering  Research  Center 

Office,  Chief  of  Engineers 

Galveston  District 

Savannah  District 

Office,  Chief  of  Engineers 

Philadelphia  District 

Coastal  Engineering  Research  Center 

San  Francisco  District 

Waterways  Experiment  Station 


Office,  Chief  of  Engineers 


Massachusetts  Institute  of  Technology 
Waterways  Experiment  Station 
Johns  Hopkins  University 
Consulting  Engineer,  Philadelphia 


FOREWORD 


i 

This  report  describes  an  analytical  approach  to  the  problem  of 
defining  changes  in  the  water  level  of  a  basin  connected  to  the  ocean 
by  a  channel  or  channels.  To  illustrate  use  of  the  analytical  approach, 
the  case  of  Indian  River,  Delaware,  is  examined  in  considerable  detail. 

The  analytical  method  was  developed  and  the  report  was  prepared  by ' ’ 
Dr.  Garbis  H.  Keulegan.  The  work  was  initiated  while  Dr.  Keulegan  was 
employed  as  Physicist  for  the  National  Bureau  of  Standards  and  was  com¬ 
pleted  while  he  was  serving  as  Staff  Consultant  to  the  Hydraulics  Divi¬ 
sion,  U.  S.  Army  Engineer  Waterways  Experiment  Station. 

The  work  reported  herein  was  sponsored  by  the  U.  S.  Army  Engineer 
Committee  on  Tidal  Hydraulics. 
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ABSTRACT 


An  analysis  is  made  of  the  changes  in  water  level  in  a  basin  con¬ 
nected  to  the  ocean  by  a  relatively  long  channel.  It  is  assumed  that  the 
banks  of  the  basin  are  vertical,  that  the  connecting  channel  is  many  times 
deeper  than  the  tidal  range,  that  there  is  no  flow  into  the  basin  from 
streams,  and  that,  therefore,  there  are  no  density  currents  present.  The 
consideration  of  the  storage  equation  shows  that  the  coefficient  of  reple¬ 
tion  is  a  characteristic  quantity  for  basins  in  general.  The  aim  of  the 
analysis  is  to  determine  the  range  of  tides  within  the  basin,  the  lag 
between  the  maximum  surface  deflections  in  the  basin  and  in  the  sea,  and 
the  maximum  mean  velocity  in  the  connecting  channel.  These  quantities  are 
determined  as  functions  of  the  coefficient  of  repletion.  Tables  are  given 
for  ready  reference  for  many  of  the  quantities. 

The  application  of  the  method  for  channels  to  flow  through  barrier 
cuts  is  also  discussed.  Finally,  the  matters  of  the  Indian  River  Inlet 
are  studied  in  detail  to  illustrate  the  practical  meaning  of  the  various 
results  of  the  main  analysis. 
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TIDAL  FLOW  IN  ENTRANCES 


WATER-LEVEL  FLUCTUATIONS  OF  BASINS  IN  COMMUNICATION  WITH  SEAS 

PART  I;  THEORETICAL  DEVELOPMENTS 

'  Introduction 

The  study  of  the  water-level  fluctuations  in  basins  communicating 
with  seas  can  be  undertaken  with  reference  to  a  number  of  conditior  s 
namely,  shape  (horizontal  projections)  of  the  basin,  slopes  of  LLv.  ■-•.Uiks, 
amount  of  inflow  from  streams,  number  of  connections  with  the  sea,  hydrau¬ 
lic  resistance  of  the  connecting  channels,  type  of  tidal  fluctuations  of 
the  seas,  relation  between  tidal  range  and  depth  of  channel,  and  presence 
or  absence  of  density  currents . 

The  movement  of  water  in  the  channel  in  the  horizontal  direction  is 
affected  in  a  very  marked  manner  by  the  shape  of  the  basin.  If  the  basin 
is  narrow,  shallow,  and  long  and  the  communication  with  the  sea  is  at  one 
end,  the  flow  of  water  from  the  sea  into  the  basin  may  be  associated  with 
advancing  waves.  In  this  case,  change  of  water  level  in  the  basin  is  not 
uniform;  hence  the  various  tidal  changes  in  the  basin  from  place  to  place 
must  be  considered.  On  the  other  hand,  if  the  basin  is  square  or  circular, 
or  of  seme  similar  form,  the  change  in  elevation  of  the  water  suiface  will 
be  nearly  the  same  for  every  point  in  the  basin. 

In  the  latter  case,  the  accumulation  of  water  in  the  basin  obeys  the 
condition  of  the  flow  of  reservoirs  and,  in  particular,  the  so-called  stor¬ 
age  law.  Hie  equation  representing  the  storage  is  also  the  differential 
equation  of  the  surface  changes.  The  iurm  of  the  equation  and  hence  the 
method  of  solving  it  are  affected  by  inflow  from  streams,  slope  of  the 
banks,  ratio  of  the  depth  of  the  connection  channel  to  the  tidal  range,  and 
way  in  which  flood  and  ebb  tides  occur. 

Hiis  report  discusses  the  simplest  set  of  conditions.  It  is  assumed 
here  that  the  walls  of  the  basin  are  vertical,  that  there  is  no  inflow  from 
streams,  that  there  are  no  density  cur1  - nts  present,  and  that  the  tidal 
fluctuations  are  given  by  a  sine  curve.  The  connecting  channel  is  assumed 


1 


to  be  prismatic,  and  1 '  .  depth  of  the  channel  is  assumed  to  be  large  with 
respect  to  the  tidal  range.  It  is  also  assumed  that  the  flow  in  the  chan¬ 
nel  is  governed  by  Manning's  formula. 

Although  this  problem  has  been  considered  by  COL  Earl  I.  Brown,  who 
has  given  a  solution,  it  was  believed  that  a  new  treatment  should  be  at¬ 
tempted  with  the  object  of  establishing  a  better  approximation.  A  novelty 
in  the  solution  presented  here  is  the  dimensionless  form  of  the  equation  of 
the  surface  changes .  The  aim  of  the  analysis  is  to  evaluate  first  the  max¬ 
imum  displacement  of  the  water  surface  in  the  basin,  and  second  the  maximum 
mean  velocity  in  the  channel  during  a  tidal  cycle. 

Derivation  of  the  Equation  of  Surface  Changes 

Let  Q  be  the  discharge  at  time  t  in  the  connecting  channel, 

V  the  mean  velocity,  Hg  the  elevation  of  the  open  sea,  and  the  ele¬ 
vation  of  the  water  surface  in  the  basin.  The  elevations  are  al 1  measured 
with  respect  to  mean  sea  level  (msl).  The  difference,  Hg  -  ,  represents 

the  fall  of  the  surface  corresponding  to  the  mean  velocity  V  (fig.  l). 


Fig.  1.  Gradient  of  water  surface  over  channel 


This  fall  can  be  broken  into  two  parts  as  follows: 


Hg  -  +  AHg 


(1) 


2 


The  first  difference  gives  the  fall  that  is  necessary  to  accelerate  the 
water  from  the  sea  entering  the  channel  from  zero  to  the  velocity  V  at 
the  entrance .  Accordingly, 

V2  ,  . 

“  m  2i  (2) 

where  m  is  a  coefficient  resulting  from  the  velocity  distribution  and  g 
is  the  acceleration  due  to  gravity.  If  the  velocity  distribution  over  the 
cross  section  is  uniform,  m  reduces  to  unity.  In  open-channel  flow  the 
exact  value  of  m  is  not  known. 

Since  the  jet  issuing  from  the  exit  end  of  the  channel  is  dissipated 
in  a  process  of  turbulent  expansion  and  the  pressures  in  the  medium  where 
the  jet  is  being  obliterated  are  hydrostatic,  the  difference  Ahg  gives 
the  fall  necessary  to  overcome  the  resistance  of  the  connecting  channel. 
Using  the  Weisbach  type  of  formula: 

<3) 

where  X  is  the  friction  coefficient,  L  is  the  length  of  the  connecting 
channel,  and  r  is  the  hydraulic  radius  of  the  channel.  The  relation  that 
exists  between  X  and  Manning's  n  will  be  considered  later. 

Combining  the  expressions  in  equations  2  and  3: 

2 

+  ^  =  (x  £+  m)g  (4) 

Introducing  the  difference,  Hg  -  ,  and  solving  for  V2  ,  yields: 

*  -  -  V  (5> 

It  is  desirable  to  express  the  >  ations  with  respect  to  msl  in  terms  of 
the  semiamplitude  of  the  tif1,l  displacements  occurring  in  the  sea.  Denot¬ 
ing  the  tidal  range  by  2H  ,  and  using  H  as  a  measure  of  the  surface 
fluctuations,  yields: 


3 


Writing  hg  =  Hg/H  ,  and  h^  =  H^/lI  , 


then 


V 


(7) 


Since  there  are  no  contributions  from  the  inflow  of  streams  and  since 
the  banks  are  vertical,  the  storage  equation  of  the  water  in  the  basin, 
assuming  that  the  tide  rises  and  falls  simultaneously  throughout  the  basin, 
is 

dlL 

AdT=sV  (8) 


where  A  is  the  surficial  area  of  the  basin,  and  a  is  the  cross- 
sectional  area  of  the  connecting  prismatic  channel.  Hence, 

=  f  V  (9) 

Denoting  the  tidal  period  by  T  permits  introduction  of  the  trans¬ 
formation 


5, 

dt 


t  _  _e 

T  2n 


(10) 


where  0  is  specific  tidal  time  in  radians,  and  equation  9  oecomes 

de  2nH  A 


(11) 


Eliminating  V  between  equations  7  and  11  and  writing 


produces  finally 


=  K^hg  - 


»  h2  >  hl 


which  is  the  differential  equation  of  the  surface  fluctuations  in  the  basin 
when  the  surface  of  the  sea  is  at  a  higher  elevation  than  the  surface  of 
the  water  in  the  basin.  When  the  condition  is  reversed,  i.e.  when  the  sur¬ 
face  of  the  sea  is  lower  than  the  surface  of  the  water  in  the  basin,  the 
corresponding  equation  is 

dh  _ 

diT  =  ~KVhl  ~  h2  9  hl  >  h2 


In  the  results  that  will  be  given  later,  the  numerical  quantity  K 
plays  a  decisive  role.  It  summarizes  the  effects  of  the  channel  and  the 
basin  dimensions,  of  the  roughness  of  the  walls,  and  of  the  period  and 
range  of  the  tidal  fluctuations  on  the  limits  of  the  water-level  changes  in 
the  basin.  Because  of  this  significance  it  appears  to  be  appropriate  to 
refer  to  K  as  the  coefficient  of  filling  or  repletion. 


aical  Values  of  the  Coefficient  of  Repletion 


It  is  desired  to  consider  some  typical  values  of  the  coefficient  K 
in  order  to  form  an  idea  of  the  variations  in  the  values  of  the  coefficient 
which  are  ordinarily  to  be  expected.  The  first  step  is  to  give  the  rela¬ 
tion  that  exists  between  the  coefficient  of  friction  X  and  Manning's  n  . 
Manning's  formula  is 


=  yl/6 

n 


where  s  ,  the  surface  slope,  is  equal  to  £Hg/L  ,  and  the  units  of  meas¬ 
urement  are  the  foot  and  the  second.  But  fran  equation  3 

V  =-»/l&""rs 
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Comparing  the  latter  two  expressions  for  the  velocity  yields 


n  V2I 

l.U86r^^ 


(15) 


which  is  the  desired  relation  connecting  X  with  n  .  Some  numerical 
values  relating  X  to  n  and  r  are  given  in  table  1. 

If  it  is  assumed  that  the  distribution  of  the  velocities  at  the  en¬ 
trance  to  the  channel  is  uniform,  m  in  equation  12  can  be  taken  as  unity. 
The  relation  giving  K  can  now  be  written  as 


AKV~H  _  T  V~2g  _  r 
a  2it  V  XL  +  r 


(12a) 


The  quantity  on  the  right  side  of  the  equation  has  been  computed  for 
various  channel  depths  and  lengths  and  for  different  values  of  n  .  Compu¬ 
tations  are  made  for  an  assumed  tidal  period  T  equal  to  12  hr.  All 
lengths  are  in  feet  and  times  in  seconds  in  the  computations .  For  the  re¬ 
sults  see  table  2.  As  an  illustration  of  the  use  of  the  table,  assume  that 
r  =  5  ft,  L  =  1000  ft,  n  =  0.03,  and  K  =  1  ft.  Then  from  table  2; 

K=  2.734  x  10k  X  | 

/  q  k  5 

so  that  if  A/a  is  given  the  succession  of  values,  10  ,  10  ,  and  10  , 

K  takes  on  the  values,  27.3,  2.73,  and  0.27,  respectively. 


Sinusoidal  Fluctuations  of  Surface  of  Sea 


Granting  that  the  fluctuations  in  level  of  the  surface  of  the  sea  can 
be  represented  as  a  pure  sine  curve,  the  most  general  solution  to  describe 
the  fluctuations  of  water-surface  level  in  the  basin  needs  to  be  effected 
with  reference  to  the  height  of  the  water  surface  in  the  basin  at  the 
instant  the  basin  is  connected  to  the  sea. 

If,  however,  the  basin  has  been  in  communication  with  the  sea  for  a 
long  time,  the  fluctuations  of  the  water  surface  in  the  basin  become  steady 
and  fluctuate  between  limits  that  do  not  vary  with  time.  Ihe  same  limits 
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are  established  no  matter  what  the  initial  depth  of  water  in  the  basin  may 
have  been.  Although  the  fluctuations  of  the  surface  of  the  water  in  the 
basin  are  periodic,  it  cannot  be  said  that  they  have  the  form  of  a  pure 
sine  curve,  because  the  frictional  resistance  of  the  connecting  channel 
varies  as  the  square  of  the  mean  velocity. 

Assume  that  the  displacements  of  the  water  surfaces  in  the  sea  and  in 
the  basin  are  given  on  a  common  axis  of  time  t  or  of  the  dimensionless 
time  parameter  6  (see  fig.  2).  The  origin  of  time  may  be  taken  as  the 


Fig.  2.  Surface  fluctuations  of  sea  and  basin 

instant  when  hg  and  h^  are  equal  and  hg  commences  to  become  greater 
than  h^  .  Then,  a's  shown  in  fig.  2,  hg  =  0  when  9=  t  .  Accordingly, 
the  oscillation  of  the  surface  of  the  sea  will  be  given  by 

hg  =  sin  (6  -  t)  ,  0  <  6  <  2n  (l6) 

A  consideration  that  will  have  a  bearing  on  the  method  of  solution  to  be 
followed  is  the  proportion  of  the  time  during  which  the  water  level  in  the 
sea  is  higher  than  the  water  level  in  the  basin.  It  will  be  assumed  that 
during  the  period  T/2  the  surface  of  the  sea  is  higher  than  the  surface 
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in  the  basin.  The  assumption  will  be  verified  later.  During  the  following 
period  the  reverse  is  true.  Accordingly,  the  determination  of  h^  will  be 
made  separately  for  the  ranges  of  values  0  <  0  <  it  and  jt  <  0  <  2it  ; 
these  ranges  will  be  referred  to  as  the  first  range  and  the  second  range, 
respectively. 

The  mathematical  task  involves  the  following  relations  for  the  first 

range : 


h2  >  hJL  ,  0  <  6  <  it 


hg  =  sin  (0  -  t) 

(16) 

,cH 

1 

V 

W 

II 

sis 

(13) 

tvf 

II 

V# 

CD 

ii 

o 

hg  =  hj  ,  0  =  it 

For  the  second  range,  the  relations  involved  are: 


and 


^  >  hg  ,  it  <  0  <  2n 


hg  =  sin  (0  -  t) 


h2  =  hi  ,  0  =  it 

h2  =  hx  ,  9  =  2it 


(16) 

(1*0 


Instead  of  determining  h^  directly,  it  is  more  convenient  to  obtain 
the  difference,  h^  -  hg  .  Thus,  putting 

z.  =  hg  -  h1  ,  0  <  0  <  «  (17) 
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equation  13  reduces  to 


dz 

d0 


After  introducing  the  value  of  h g  from  equation  16. 
cal  problem  for  the  first  range  becomes  that  of  determining 
following  relations : 


z  >  0  ,  0  <  0  <  it 


t  ' 

=  -K  -y  z  +  cos  0  cos  t  +  sin  0  sin  t 


z  =  0  ,  0=0 


z  =  0  j  0  =  « 


For  the  second  range,  putting 


z  =  h^  -  hg  ,  it  <  0  <  2it 


equation  14  becomes 


S=-Kyi- 


d0 


dhg 

dlT 


Introducing  transformation 


0  =  «  +  P 


and  hence  the  relation 


hg  =  -sin  (0  -  t) 

the  mathematical  problem  for  the  second  range  involves  the 
conditions 


o 

s 


the  mathemafci- 
z  from  the 


(18) 


(19) 


(20) 


solution  of  the 


z>0,0<P<Jt 


•ts-  -  -K  V z  +  cos  p  cos  t  +  sin  P  sin  t 
dp 


z  =  0  ,  P  =  0 


■21) 


z  =  0  j  P  =  n 

Comparison  of  the  systems  of  relations  in  equations  18  and  21  indi¬ 
cates  that  it  should  be  sufficient  to  obtain  the  solution  of  z  for  the 
first  range.  The  behavior  of  z  for  the  second  range  is  readily  deduced. 
The  same  comparison  shows  also  that  the  portion  of  the  time  during  which 
the  surface  of  the  sea  is  at  a  higher  level  than  the  surface  in  the  basin 
is  T/2  .  Thus  the  original  assumption  is  confirmed. 

Now  since  the  solution  of  equation  18  is  of  the  form 

z  =  f(6)  (22) 


and  since 


h1  =  h2-z,O<0<Jt 
the  value  of  h^  for  this  range  is 

hx  -  sin  (0  -  t)  -  f(0)  (23) 

Again,  since  the  solution  of  equation  21  is  of  the  form 

z  =  f(p)  (24) 


and  since 


hl=z  +  h2,°<P<n 
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hx  =  sin  0  -  t)  +  f(p)  ,  0  =  it  +  P  (25) 

The  geometrical  interpretation  of  the  above  discussion  is  readily  in¬ 
ferred.  The  curve  of  h^  for  the  first  range  is  drawn.  The  curve  is  re¬ 
flected  and  is  moved  along  the  axis  of  0  by  an  amount  jt  .  The  curve 
thus  displaced  is  the  curve  of  h^  for  the  second  range. 

Fourier  Expansion  of  Vsin  0 

The  determination  of  a  form  of  z  that  will  satisfy  the  differential 
equation,  equation  18,  can  be  effected  in  various  ways.  The  solution  can 
be  in  terms  of  polynomials  of  0  or  in  terms  of  the  circular  functions  of 
9  .  Since  periodic  changes  are  involved,  it  is  preferable  to  obtain  the 
solution  in  circular  functions.  Select  as  a  possible  expression  the  series 


z  =  51  An  sin  n0  +  ^  B^fcos  nd  -  cos  (n  +  2)0] 


n' 


1  1 

n  =  1,  3>  5»  ...»  2m  +  1 


(26) 


since  z  vanishes  for  6  =  0  and  9  =  n  .  Here,  A  and  B  are 
constants . 

It  is  obvious  that  the  first  term  in  the  series,  i.e.  the  quantity 
A^  sin  0  ,  is  the  predominant  term.  In  equation  18  there  is  the  term  ~fz  , 
and  it  will  be  necessary  in  the  course  of  the  analysis,  as  will  be  seen 
later,  to  have  the  Fourier  expansion  of  "V  sin  9  .  Now,  by  the  rule  of  the 

Fourier  expansion,  if  f ,  ( 0)  is  single -valued,  finite,  and  continuous  be- 

i  2 
tween  0=0  and  0  =  n  ,  it  can  be  developed  into  a  series  of  the  form 

^(0)  =  a^^  sin  0  +  &2  sin  20  +  a^  sin  30 
where  the  coefficients  have  the  values 
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I 


m 


2  r51 

=  —  /  f^(0)  sin  ra©  d0 


In  conformity  with  the  rule  and  by  inspection,  the  following  equation  can 
be  written: 


l/sin  0  =  sin  0  +  sin  36  +  fy.  sin  5 0 


(27) 


where 


v-Sfv= 

N3  -  C  V 

“5  -  C  VH 


sin  6  sin  6  d 0 


sin  0  sin  30  d0 


sin  0  sin  50  d0 


(28) 


The  N's  can  be  determined  numerically  by  replacing  the  process  of 
integration  by  a  process  of  summation.  For  example: 


~^s in  6  sin  0  AS  ,  0  <  0  <  n 


The  summations  are  made  by  putting  the  interval  AS  equal  to  0.03^91 
radians.  The  computations  give 


Nx  =  1.1107 


n3  =  0.1580 


Nc  =  0.0711 

5 


(29) 


The  degree  of  convergence  of  the  series  in  equation  27  will  be  exam¬ 
ined  best  when  the  differences 


A^  =  ~\fsin  0  -  sin  0 
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=  "ysin  Q  -  Nx  sin  6  -  N3  sin  36 

^5  =  Vsin  9  ~  N!  sin  0  “  N3  sin  30  -  N,_  sin  5© 

are  computed  and  plotted  (see  fig.  3)*  The  convergence  appears  to  be  sat¬ 
isfactory.  The  larger  deviations  are  confined  to  the  end  regions.  The 
effect  of  adding  sine  terms  of  the  higher  frequencies  is  to  decrease  the 
interval  length  of  the  end  regions  where  the  maximum  deviations  occu^ .  In 


0,  DEGREES 

Fig.  3.  Residuals  in  the  successive  approximations  of  the 
Fourier  expansion  of  ysin  e 


addition,  the  magnitudes  of  the  deviations  are  decreased.  It  appears  also 
that  for  the  ordinary  approximations,  when  very  high  accuracy  is  not  needed, 
it  is  sufficient  to  adopt: 

~\J sin  0  =  sin  0  +  sin  30  (30) 

Conforming  to  this  selection,  the  terms  of  higher  frequencies  in  the  ex¬ 
pression  for  z  ,  equation  26,  could  be  dispensed  with. 

Method  of  Approximate  Solutions 

As  was  mentioned  above,  it  will  be  sufficient  to  determine  z  for 
the  first  range  only,  i.e.  for  the  interval  0  <  0  <  it  .  The  differential 
equation  to  be  solved  is 

cLz  r™ 

rj'g  =  -K  y  z  +  cos  0  cos  r  +  sin  0  sin  x  (l8) 

The  solution  is  taken  to  be  of  the  form 


z  =  a^  sin  0  +  a^b^fcos  0  -  cos  3©)  +  sin  30 


(31) 


an  expression  which  vanishes  for  the  points,  0=0  and  0  =  r  .  Since 
terms  of  higher  frequency  are  omitted,  the  solution  represents  an  approxi¬ 
mation.  The  unknown  quantities  are  the  coefficients,  a^  ,  a^  ,  b^  ,  and 
the  phase  angle  x  . 

Factoring  out  the  first  term  on  the  right  side  of  equat:  on  31  yields 


z  = 


*1 


sin 


#C: 


1  +  b. 


cos  0  -  cos  30 
sin  0 


+  a. 


sin  30 
sin  0 
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(32) 


Since  the  terms  in  parentheses  remain  finite  when  0  is  varied  from  0  to 
it  ,  it  is  admissible  to  take  the  square  root  of  both  sides.  Hence, 


z 1'2  =  a^/2  ~\j  sin  0  ^1  + 


b3  cos  0  -  cos  30  , 
2  sin  0 


^3  sin  30^ 
2  sin  0  J 


(33) 
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It  is  conceived  that  and  are  fractions;  therefore,  in  writing  the 
square-root  expression  for  the  terras  in  the  parentheses  of  equation  32,  the 
terms  containing  the  squares  of  a^  and  bg  and  their  product  are  ne¬ 
glected.  Introducing  the  Fourier  expression  of  equation  30  yields 


s1/2  =  af  (% 


sin  0  +  sin 


36)  ( 


1  + 


h 


cos  0  - 


sin 


cos 

V 


36 


+ 


1  sin  30\ 

2  3  sin  e  J 


(34) 


Removing  the  first  set  of  parentheses  and  ignoring  the  terms  multiplied  by 
N^b^  and  N.^  ,  since  these  are  small  fractions,  yields 


z1/2  =  a^/2  sin  0  +  sin  30  +  -g1 ' (cos  0  -  cos  3e) 


a  ^ 

+  -g —  sin  30 


(35) 


Differentiating  equation  31  yields 


~  =  a.^  cos  0  +  a^b^(-sin  9+3  sin  39)  +  Sa^a^  cos  30  (36} 


Substituting  the  expressions  from  equations  35  and  36  in  equation  18, 
the  latter  equation  reduces  to 

^ a^/2  N^K  -  a.^  -  sin  sin  0  +  ^a*/2  b^  +  a1  -  cos  cos  0 
+  j^/2  K  (n3  +  —  aj  +  3a1b3j  sin  3 0 

+  (_aJ-/2  ^  b3  +  3^)005  30  =  0  (37) 


This  equation  must  hold  for  any  value  of  0  .  Hence  it  must  hold 
when  multiplied  by  sin  0  d 0  and  integrated  between  the  limits  0  and 
Jt  .  The  same  is  true  when  it  is  multiplied  by  cos  0  d0  ,  cos  30  d0  ,  or 
sin  30  d0  and  integrated  between  the  same  limits.  Carrying  out  all  the 
steps,  it  is  found  that 
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(38) 
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r~ 


N  I  N^K2  \ 

a3=-2^(^?71^)  (47) 

a/2 

*3  “  “3  <«> 

With  the  values  of  and  b^  thus  known,  the  second  approximate 

values  of  a^  and  t  can  be  obtained  by  reverting  once  more  to  equations 
38  and  39.  Writing 


it  is  found  that: 


®i  =  ®n  +  5ai 


T  =  +  6T 


«K2  v 


8i-  -  5  + 


Bearing  in  mind  that  the  fluctuations  of  the  water  surface  in  the 
basin  are  given  by  h^  =  z  +  h^  and  that  z  is  given  by  equation  31j  the 
quantities  a_^  ,  a^  ,  b^  ,  and  t  are  the  parameters  that  determine  the 
form  of  the  fluctuations  in  the  basin  as  a  function  of  time.  These  con¬ 
stants  depend  individually  on  the  coefficient  of  repletion  K  .  Their 
values,  decermined  according  to  the  scheme  of  relations  discussed  above, 
are  given  in  table  3* 

In  the  method  used,  the  results  are  obtained  within  the  solution  to  a 
second  approximation  of  the  differential  equation,  equation  18.  The  method 
of  analysis  is  such  that  one  can  go  to  higher  approximations.  This  possi¬ 
bility,  however,  has  only  a  theoretical  significance,  since  the  computa¬ 
tions  to  be  made  are  very  long  and  should  be  avoided. 
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Surface  Displacement  Curves  for  the  Case  when  K  is  Equal  to  Unity 

By  way  of  illustration,  the  surface  curves  for  the  case  in  which  the 
coefficient  of  repletion,  K  of  equation  12,  is  equal  to  unity  will  be  de¬ 
termined.  Tab3e  3  shows  that  in  this  case  the  surface  curve  parameters  are: 

a1  =  0.5^51  cos  r  =  0.5178 

a^  =  -O.OI65  sin  t  =  0.8555 

b3  =  -0.0664 


Substituting  these  values  in  equations  16  and  31,  the  quantities  hg 
and  z  are  determined.  Since  z  =  h^  -  hg  ,  h^  is  computed  by  taking  the 
sum  of  z  and  hg  .  The  results  are  shown  in  fig.  4.  In  preparing  the 


Fig.  4.  Surface  fluctuations  ,of  sea  and  basin  waters  for  the  case  K  =  1 

f 

plots,  use  was  made  of  the  fact  that  the  solutions  for  the  range  0  <  6  <  « 
can  be  extended  to  the  range  n  <  9  <  2n  .  This  matter  has  been  discussed 
previously. 
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The  solutions  given  are  only  approximate,  and  it  would  he  instructive 
to  find  how  closely  the  differential  equation,  equation  18,  is  satisfied. 
Because  of  the  approximations,  if  the  derived  solution  is  reintroduced  into 
the  differential  equation,  a  remainder  AR  will  be  left;  that  is, 

^  =  if +  K  V*  -  w  (53) 

The  smaller  the  difference  AR  ,  the  closer  the  approximation  will  be 
to  the  exact  solution.  In  this  case  an  idea  can  be  obtained  concerning  the 
sufficiency  of  the  approximation  by  comparing  the  remainder  AR  with 
dh2/de  .  For  this  comparison  the  quantities  AR  and  dhg/d©  are  plotted 
against  6  in  fig.  5.  It  can  be  seen  that  the  remainder  is  small,  and 
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that  a  good  solution  has  been  obtained  in  the  case  under  consideration. 

Range  of  Tides  in  the  Basin  and  Lag  of  the  Maxima 

The  range  of  the  tides  in  the  basin  is  twice  the  maximum  displacement 
of  the  water  surface  in  the  basin  measured  from  msl.  The  maximum  and  the 
minimum  displacements  correspond  to  the  zeroes  of  the  difference  z  (i.e. 
h5  =  h^),  since  for  these  points  the  rate  dh^/de  vanished  (see  equation 
13).  As  z  vanishes  at  the  points  0=0  and  0  =  it  ,  it  is  sufficient 
to  consider  the  value  of  h^  at  0  =  it  ,  where  h^  is  maximum.  Let  this 
value  be  h1m  .  Bearing  in  mind  the  dimensionless  character  of  the  quan¬ 
tity,  it  is  seen  that  h^  gives  the  ratio  of  the  semirange  of  tides  in 
the  basin  to  the  semirange  of  tides  in  the  sea.  Since  at  0  =  it  ,  h^ 
equals  Ip,  ,  and  the  value  of  hg  at  0  =  it  is  sin  r  (see  equation  16), 
the  ratio  of  the  semirange  of  tide  in  the  basin  to  the  semirange  of  tide  in 
the  sea  is 


h^  =  sin  t  (54) 

The  ratio  of  the  range  of  tide  in  the  basin  to  the  range  in  the  sea  is  also 
sin  t  .  The  values  of  sin  t  as  a  function  of  K  are  shown  in  table  4. 
Thus  the  tidal  range  of  the  water  in  the  basin  can  be  read  directly  from 
the  table,  once  the  coefficient  of  repletion  K  is  known  for  the  particu¬ 
lar  basin  (see  also  fig.  6). 

The  next  question  to  consider  is  the  lag  between  the  maximum  dis¬ 
placement  of  the  water  surface  of  the  sea  and  the  maximum  displacement  of 
the  water  surface  in  the  basin.  In  order  of  time,  the  former  precedes  the 
latter.  Let  the  lag,  expressed  in  radians,  be  denoted  by  a  .  The  maximum 

displacement  of  the  surface  of  the  sea  occurs  at  0  and  has  the  value, 

m 

from  equation  16, 


0 


m 


T 


2 


T 


The  maximum  displacement  of  the  surface  of  the  water  in  the  basin  occurs 
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coefficient  of  repletion 


when  0  -  «  .  Hence  the  lag  is : 

a  =  n  -  6 
m 

or 

<*  =  §-t  (55) 

The  question  of  sin  t  was  considered  in  the  preceding  paragraph. 


Tidal  Prisms  and  Maximum  Mean  Velocity  in  the  Communicating  Channel 


The  volume  of  water  in  the  basin  included  between  the  two  horizontal 
planes,  one  giving  the  highest  elevation  of  the  water  surface  during  a 
tidal  cycle  and  the  other  giving  the  lowest  elevation  of  the  surface,  is 
referred  to  as  the  tidal  prism.  Let  the  volume  of  the  prism  be  ft  .  If 
is  the  maximum  rate  of  discharge  through  the  connecting  channel  during 
a  half  tidal  cycle,  the  volume  of  the  prism,  the  maximum  rate  of  discharge, 
and  the  period  of  the  tides  can  be  connected  by  an  expression: 


3» 

An 


=  c 


(56) 


where  C  is  a  dimensionless  number. 

The  value  of  C  is  close  to  unity,  and  its  exact  value  depends  on 
the  coefficient  of  repletion  K  .  This  dependence  will  be  determined  next. 
With  the  maximum  mean  velocity  being  denoted  by  Vffi  , 


Also,  from  the  condition  of  continuity. 


Since  z  =  hg  -  h^  ,  it  can  be  seen  from  equation  7  that 


V~  Vz 
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and 


v  ~  (V^L 

m  ’  m 


the  suffix  m  indicating  that  maximum  values  are  taken.  Hence, 


or,  since  2ndt  --  Td0  , 


(57) 


Let  8^  be  the  value  of  G  which  renders  z  a  maximum.  At  this 
point  r*  is  also  a  maximum.  Accordingly,  dz/d0  =  0  ,  and  from  equa¬ 
tion  36: 


cos  +  b^(-sin  6^+3  sin  36^)  +  3a^  cos  36^  =  0 
It  can  be  shown  that  the  smallest  root  of  this  equation  is: 


e 


(58) 


where 

4 

6  “l-9a3*3 

and 


cos  0^  =  -e 
cos  3&j_  =  +3e 
sin  6^  =  1  ,  sin  3^  =  -1 
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Prom  equation  35  the  maximum  value  of  v*  is 

(Vz)m  =  ^  sin  ei  +  N3  sin  3^ 

b-N.  Na  I 

+  (cos  91  -  cos  3^)  +  — sin  36-J 

Introducing  the  value  of  e1  from  the  preceding  page,  this  reduces 

(V^)m  =  aj/2  (n±  -  N3  -  2Nxb3e  -  (59) 

Again,  from  equation  35 ,  effecting  an  integration: 


z^2de  = 


-  2^/a  (B1  +  5  n3  +  I  Nias) 


Substituting  these  expressions,  equations  59  and  60,  in  equation  57, 
and  then  making  use  of  the  fact  that  N3  ,  a3  ,  b3  ,  and  r  are  all  small 
quantities,  yields 

Jtft  3  Nx  3  3  3 

The  right  member  of  the  equation  is  the  expression  for  C  appearing  in 
equation  56;  that  is, 


C  1  '  3^'  2b3e  "  3  a3  (61' 

It  is  obvious  that  C  depends  on  K  ,  since  a3  ,  b3  ,  and  e  de¬ 
pend  on  K  .  Values  of  C  computed  using  the  latter  expression  are  given 
in  table  4.  It  is  seen  that  as  K  is  increased  from  0.1  to  100  the  value 
of  C  changes  from  0.8106  to  1.0000  (see  also  fig.  6). 

The  formula  for  the  tidal  prism  is  the  means  by  which  the  maximum 
velocity  in  the  connecting  channel  can  now  be  evaluated.  By  the  definition 
of  the  tidal  prism, 
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where  is  the  ratio  of  the  maximum  displacement  of  the  water  surface 

in  the  basin  to  the  maximum  displacement  of  the  surface  of  the  sea.  Hence, 
using  equation  54 


ft  =  2AH  sin  r 


As  before, 


Introducing  these  into  the  formula  for  the  tidal  prism,  equation  56, 

gives 

V  =  2nC  -  J  sin  r  (62) 

m  a  T  '  ' 

This  is  the  relation  which  connects  the  maximum  mean  velocity  in  the 
connecting  channel  with  the  range  of  tides  in  the  sea,  2H  .  In  the  ex¬ 
pression  K  is  measured  in  feet  and  T  in  seconds. 


Prismatic  Equivalence  of  Irregular  Channels 
and  Multiplicity  of  Channels 

The  determination  of  the  coefficient  of  repletion  K  of  the  connect¬ 
ing  channel  has  been  carried  out,  as  discussed  in  the  initial  sections  of 
this  report,  on  the  basis  that  the  channel  is  prismatic.  An  actual  channel 
may  not  meet  this  condition,  the  depth  and  width  varying  from  point  to 
point.  For  the  purpose  of  evaluating  K  ,  the  actual  irregular  channel  can 
be  replaced  by  a  regular  one  having  the  same  conductance  as  the  irregular 
channel.  In  finding  the  equivalent  channel  it  may  be  well  to  assume  that 
the  channel  length  L  and  the  roughness  coefficient  n  remain  the  same; 
also  that  the  depth  of  the  equivalent  channel  is  equal  to  the  average  depth 

r  of  the  irregular  channel.  Thus  it  remains  to  find  the  constant  cross 
m 

section  a  of  the  equivalent  channel.  If  the  end  cross  sections  of  the 
s 

actual  channel  are  not  of  the  same  size,  the  activity  of  the  channel 
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differs  for  the  two  directions  of  flow,  one  direction  being  the  reverse  of 
the  othi‘r.  If  then  agl  is  the  equivalent  cross  section  for  one  direction 
and  ag2  for  the  opposite  direction,  it  should  suffice  to  take  the  mean 
(asl  +  as2^  as  the  cross  secbion  of  the  equivalent  channel.  If  the 
length  of  the  connecting  channel  is  large  in  comparison  with  the  depth,  the 
error  arising  from  this  difficulty  will  be  reduced. 

Let  r  and  a  be  the  depth  and  the  cross-sectional  area  of  the 

X  X 

irregular  channel  at  the  point  x  .  Let  a  be  the  mean  value  of  a  . 

m  x  ’ 

averaged  along  the  entire  length  of  the  channel.  Let  a^  and  a2  be  the 
end  cross  sections  of  the  channel,  a^  being  on  the  end  toward  the  sea. 
Assume  that  the  flow  Q  is  from  the  sea  toward  the  basin  (see  fig.  7). 


SEA 


^77777777777777777777777777777777777777777777777^ 


I 


BASIN 


Fig.  7.  An  irregular  connecting  channel 


The  surface  fall  may  be  broken  into  two  parts,  and  AHg  ,  having 

the  values 
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2  X  L  n2 

“  2g  2  2g  I  2 
B  r  a  6  a 

m  s 

In  these  expressions  X^  is  the  coefficient  of  friction  correspond¬ 
ing  to  the  hydraulic  radius  r^  and  n  ,  and  X^  is  the  coefficient  of 
friction  corresponding  to  the  hydraulic  radius  r^  and  n  .  Since  chan¬ 
nels  are  very  wide,  the  average  depth  of  a  cross  section  can  be  taken  equal 

to  the  hydraulic  radius  of  the  section.  In  channels  having  equal  conduct- 

2 

ances,  the  quantity  2gAH/Q  is  the  same  for  the  same  value  of  Q  . 
Applying  this  to  the  above  gives 


1_ 

2 

a 


If  it  is  assumed  that  along  the  channel  the  variation  in  depth  is  not 

very  great,  the  above  expression  can  be  simplified.  Putting  rx  =  rffi  +  6r  , 

it  can  be  inferred  from  equation  15,  since  6r  is  small  with  respect  to 

r  ,  that 
m 


Substituting  this  in  the  above  equation,  multiplying  the  resulting 

2  XjjjL 

equation  by  affl  ,  and  then  dividing  by  -j—  +  1  yields 

m 


which  is  the  formula  for  determining  the  cross  section  of  the  equivalent 
channel  when  the  flow  is  frcm  the  sea  to  the  basin.  If  a^  is  replaced  by 
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i 
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* 
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a1  in  equation  63,  the  result  gives  the  cross  section  of  the  equivalent 


channel  for  flow  in  the  other  direction. 

When  there  is  more  than  one  connecting  channel,  it  is  a  simple  matter 
to  determine  the  equivalent  single  channel.  In  fact,  if  ,  Kg  ,  , 


...  are  the  coefficients  of  repletion  of  the  individual  channels,  the  co¬ 
efficient  of  repletion  of  the  single  equivalent  channel  is 


K=K1+K2+K3  +  *' 


m 


This  fact  will  be  readily  understood  if  the  derivation  of  equation  13 
is  borne  in  mind. 


Empirical  Values  of  the  Tidal  Prism  Formula 


Keulegan  and  Hall  have  given  an  empirical  determination  of  the  con- 
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stant  C  in  the  tidal  prism  formula.  Based  on  the  observed  values  of  the 
actual  discharges  for  the  inlets  at  Nantucket,  Manasquan,  Beaufort,  and 
Baker's  Haulover,  the  average  value  of  C  is  found  to  be  0.86. 

Considering  further  the  tidal  ranges  in  the  basins  involving  these 
inlets  and  the  ranges  of  the  seas  outside  during  the  period  of  the  dis¬ 
charge  observations,  it  will  be  possible  to  obtain  theoretical  values  of 
C  according  to  the  methods  discussed  previously.  In  table  4  the  quanti¬ 


ties  h1m  (i.e.  sin  t)  and  C  are  given.  Thus  it  will  suffice  to  take  the 


ratios  of  the  tidal  ranges  in  the  basins  and  the  seas  and  deduce  therefrom 
by  referring  tc  table  4  the  corresponding  values  of  C  .  This  is  shown 
below  for  the  inlets  under  consideration. 


Inlet 

Inside  Outside 

h7 

1m 

C 

Nantucket 

2.0 

2.5 

0.800 

0.841 

Manasquan 

2.5 

3.7 

0.676 

0.825 

Beaufort 

3.5 

4.0 

O.878 

0.862 

Baker's  Haulover 

3-0 

3-3 

0.910 

0.872 

For  these  four  inlets  the  average  theoretical  value  of  C  is  O.85, 
and  this  compares  well  with  the  empirically  determined  value. 
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Despite  the  irregular  and  rugged  shape  of  the  basin  areas  of  the 
respective  basins,  the  above  agreement  would  indicate  that  probably  in  the 
basins  considered  the  total  variations  are  nearly  constant  over  the  entire 
areas  of  the  basins. 

Comparisons  with  COL  Brown’s  Analysis 

The  original  idea  that  tht  flow  through  an  inlet  is  essentially  hy¬ 
draulic  stemmed  from  COL  Brown.  The  problem  was  approached  by  assuming 
that  if  the  surface  fluctuations  in  the  sea  are  sinusoidal,  then  the  fluc¬ 
tuations  of  the  lagoon  also  are  sinusoidal.  Since  the  resistance  of  the 
inlet  channel  is  best  described  by  quadratic  laws,  the  fluctuations  in  the 
lagoon  fail  to  conform  to  a  pure  sinusoidal  variation.  The  developments  on 
the  preceding  pages  were  meant  to  improve  the  rigor  of  analysis.  There¬ 
fore,  it  may  be  worthwhile  to  make  a  comparison  between  the  method  reported 
herein  and  those  of  COL  Brown's  analysis  for  a  hypothetical  basin.  In  the 
review  of  this  report  Mr.  Wicker  made  this  comparison,  and  the  results  are 
reproduced  below  from  the  commentary  kindly  offered. 

Assume :  A  basin  with  A  =  hOO  X  10^  ft2 
A  channel  with  r  ®  lp  f t 

a  =  10,000  ft2 
L  »  5,000  ft 
n  =  0.04 
An  ocean  tide  2H  =  4  ft 
H  =  2  ft 

Then,  according  to  the  procedure  reported  herein,  from  table  2 

— ^  X  io-i+  =  2  .o4o 
a 

a no  thus 

K  =  0.36 

From  table  4  for  K  =  O.36  (see  also  fig.  6) 

sin  x  =  0.40  and  C  =  0.8l4 
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Then  from  equation  62 


Vm  =  3-78  ft/sec 

i 

j 

According  to  the  procedure  of  COL  Brown,  * 


0,  =  Ah 

=  17-0kk a  Vm  “  17-0l,1,s  ^  rV3VS  ‘\/?  -  n2 


where  notation  is  the  same  as  used  herein  except^  that 
Qj.  -  total  inflow  (or  outflow) 

H  =  range  of  tide  in  sea;  i.e.  2H  in  ■ 

the  notation  of  this  report  j 

h  =  range  of  tide  in  basin;  I 

i.e.  2H.  in  the  notation  of  this  \ 

report  j 

» 

Assuming  that  h  =  1.8  ,  j 

=  720  X  106  ft3  j 

Also,  from  the  second  expression  for  in  the  above, 

Qt  =  730  X  106  ft3 

Assuming  that  h  =  1.9  > 

Qj.  =  760  X  lo6  ft3 

Qj.  =  723  X  106  ft3 
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V  =  3.52  ft/sec 
m  ' 

According  to  the  procedure  of  COL  Brown,  assuming  that  h  =  2  , 

(^  =  800  X  106  ft3 
Q.  =  329  X  106  ft3 

Assuming  that  h  =  1  , 

Qj.  =  400  X  106  ft3 
^  =  337  X  106  ft3 

Assuming  that  h  =  0.8  , 

Qj.  =  320  x  106  ft3 
^  =  338  X  106  ft5 

Evidently  h  will  he  approximately  0.8.  Taking  this  value, 

Vm  =  3.96  ft/sec 

which  compares  with  the  value  of  3*52  from  the  procedure  reported  herein. 
COL  Brown’s  procedure  again  overestimates  as  compared  with  the  results  of 
the  method  herein.  The  difference  is  about  the  same  as  previously  computed. 

Application  to  Barrier  Cuts 

The  procedure  of  evaluations  relating  to  inlets,  as  discussed  previ¬ 
ously,  is  equally  applicable  to  barrier  cuts  separating  two  bays  or  the  two 
parts  of  a  bay,  provided  the  coefficient  of  repletion  K  is  determined 
anew. 
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Let  Q  be  the  discharge  at  the  cut,  V  the  mean  velocity,  the 
elevation  of  water  surface  on  one  side  of  the  cut,  and  the  elevation 

in  the  protected  bay.  As  before,  the  elevations  are  measured  with  respect 
to  msl.  The  flow  through  the  cut  will  be  thought  to  be  analogous  to  flow 
through  an  orifice  connecting  two  reservoirs.  Then,  let  be  the  coef¬ 
ficient  of  discharge  and  a  the  mean  cross  section  of  the  current  at  the 
cut .  Hence 


Q  =  Cda  -  1^) 


(65) 


Using  the  half  tidal  range  H  in  the  outer  bay  as  a  measure  of  surface  de¬ 
flections, 


V  = 


(66) 


Using  the  surface  area  of  the  protected  bay  A  ,  the  storage  equation  is 


^l__T_av 

d©  2*H  A  V 


(ID 


Eliminating  V  between  equations  11  and  66  yields 

/ - 

d6  ^  "V  h 2  ”  hl  *  **2  >  hl 

d^  i - 

d©  “  -KL  y  hl  "  h2  ’  hl  >  h2 


(13) 

(14) 


where  the  coefficient  of  repletion  now  has  the  new  value 


jr  "JIM: 

T.  2nH 


_T  a  r 
A  Cd 


(67) 


Once  is  obtained,  the  values  of  C  and  sin  t  are  read  from  table  4 
or  from  fig.  6.  Using  these  determinations,  the  maximum  water  elevation  in 
the  protected  part  of  the  bay,  the  maximum  mean  current  through  the  cut, 
and  the  lag  in  the  oscillations  in  the  protected  part  are  established 
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I 


from  equations  5k,  62,  and  55>  respectively. 

Similar  to  the  case  of  a  canal,  assuming  that  the  surface  oscilla¬ 
tions  on  the  two  sides  of  the  cut  for  points  sufficiently  removed  from  the 
cut  are  nearly  sinusoidal,  a  would  be  the  cross  section  of  the  cut  when 
the  waters  on  the  two  sides  of  the  cut  are  at  msl.  As  regards  the  coef¬ 
ficient  of  discharge  ,  ignoring  the  effects  of  the  geometry  of  the  cut 
and  the  roughness  of  the  passage,  this  may  be  put  equal  to  unity. 


3^ 
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PART  II:  APPLICATION  TO  INDIAN  RIVER  BAY 

Topographic  Conditions  of  Indian  River  Inlet 

At  the  suggestion  of  the  Committee  on  Tidal  Hydraulics,  the  case  of 
Indian  River  Bay,  Del.,  was  examined  in  detail  to  illustrate  the  applica¬ 
tion  of  the  formulas  of  the  present  analysis. 

Indian  River  Bay  is  a  lagoon  about  2  miles  wide  from  north  to  south 
and  6  miles  long,  the  surface  area  being  15  square  miles. ^  The  main  trib¬ 
utary  is  the  Indian  River,  which  flows  into  the  bay  from  the  west.  For 
bout  3-1/2  miles  above  its  mouth  the  river  is  broad  and  relatively  shal- 
1.  .  having  a  depth  of  about  7  ft  at  local  mean  low  water  (mlw).  From 
this  point  to  the  head  of  navigation,  a  distance  of  about  8  miles,  it  is 
narrow  and  sluggish  with  a  meandering  course,  bordered  by  mud  flats  and 
swamps.  Depth  of  water  in  the  bay  averages  about  3  ft  at  local  mlw.  The 
bay  is  connected  with  the  ocean  through  the  Indian  River  Inlet. 

This  bay  is  connected  with  Rehoboth  Bay  through  two  waterways  known 
as  "The  Ditches,"  see  fig.  8.  These  passages  had  a  controlling  depth  of 
only  2  or  3  Ft  below  the  local  mlw  prior  to  1956. 

Rehoboth  Bay  is  a  shallow-bordered  lagoon  separated  from  the  ocean 
by  a  barrier  beach.  It  is  3  ailes  wide  from  east  to  west  and  4  miles  long, 
with  a  surface  area  of  approximately  13  square  miles.  The  controlling 
depth  is  3  ft  below  mlw.  Two  tributary  streams,  Herring  Creek  and  Love 
Creek,  enter  the  bay  from  the  west.  Lewes  and  Rehoboth  Canal  is  an  inland 
waterway  from  Rehoboth  Bay  to  Delaware  Bay.  From  the  jettied  entrance  on 
the  north  shore  the  inland  waterway  follows  a  northerly  course  through 
marshland,  entering  fastland  in  the  vicinity  of  Rehoboth,  then  in  a  north¬ 
westerly  direction  it  again  passes  through  marshland  to  an  inlet  from  Del¬ 
aware  Bay,  and  thence  to  a  connection  with  Broadkill  River.  Prior  to  1937 
the  total  length  was  11.8  miles.  The  project  depth  is  6  ft  below  local 
mlw  and  bottom  width  varies  from  50  to  100  ft.  It  must  be  kept  in  mind 
that  the  connection  of  Lewes  and  Rehoboth  Canal  with  Delaware  Bay  via  Roo¬ 
sevelt  Inlet  was  not  in  existence  in  1928-1929.  The  length  given  here  is 
for  a  connection  with  Delaware  Bay  via  Broadkill  River. 
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VICINITY  MAP 
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LOCATION  MAP 


Fig.  8.  Location  and  vicinity  maps  of  Indian  River  Bay  and  Rehoboth  Bay 
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Indian  River  Inlet 


5 

Indian  River  Inlet  has  an  interesting  history.  Prior  to  the  mod¬ 
ern  improvements  it  was  an  unstable  passage,  opening  and  closing  alter¬ 
nately.  Under  such  conditions  an  inlet  is  likely  to  move  from  one  place 
to  another,  and  evidence  suggests  that  during  the  period  from  1847  to  1910 
this  inlet  migrated  northward  a  distance  of  2  miles.  About  1910  the  nat¬ 
ural  inlet  had  shoaled  to  an  extent  that  a  closing  was  impending.  In 
fact,  about  the  year  1925  the  closing  was  complete  and  little  hope  could 
be  entertained  for  the  opening  under  natural  causes.  In  the  summer  of 
1928  the  Indian  River  Commission  undertook  to  open  the  inlet  by  making  a 
cut  through  the  barrier  reef  at  a  point  almost  opposite  the  center  of 
Indian  River  Bay  and  1  mile  south  of  the  final  position  of  the  natural 
inlet.  A  channel  60  ft  wide  and  4  ft  deep  at  mlw  of  the  ocean  was  exca¬ 
vated  from  the  bay,  beginning  where  the  water  was  4  ft  deep  on  the  same 
datum,  extending  eastward  through  the  marsh,  and  then  partly  through  the 
barrier  reef  to  within  115  ft  of  the  high-water  mark  on  the  ocean  beach. 
Actual  work  began  on  23  June  1928  and  was  finished  in  October.  From  the 
eastern  end  of  the  dredged  cut  an  attempt  was  made  to  complete  the  open¬ 
ing  to  the  ocean  by  first  removing  a  part  of  the  sand  with  scrapers  and 
then  blasting  out  a  channel.  The  final  attempt  in  November  1928  was  not 
successful,  and  as  the  funds  were  exhausted,  work  was  suspended.  About  a 
month  after  the  work  had  stopped  the  locality  regained  its  former  appear¬ 
ance  except  for  the  opening  through  the  dunes  and  the  cut  in  the  marsh. 

The  trench  made  by  the  explosions  was  obliterated.  The  ocean  end  had.  been 
filled  to  above  low  water  for  a  distance  of  125  ft.  At  the  end  of  April 
1929,  locally  interested  people  made  a  small  cut,  and  the  impounded  water 
eroded  an  opening  through  the  remainder  of  the  barrier  beach.  The  cross 
section  of  this  opening  was  estimated  to  be  100  ft  wide  and  4-1/2  ft  deep 
in  June  1929,  but  late  that  summer  it  was  evident  that  the  inlet  was  shoal¬ 
ing  rapidly  and  again  was  about  to  close  entirely.  In  November  1929, 
dredging  of  a  channel  60  ft  wide  at  the  bottom  and  8  ft  deep  at  mlw  was 
commenced,  beginning  in  the  marsh  about  3700  ft  from  the  ocean.  The  ex¬ 
cavation  followed  the  line  of  the  previous  cuts,  and  by  January  1930  it 
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had  reached  the  outer  bar.  The  work  provided  an  inlet  that  was  con¬ 
sidered  to  be  comparable  in  size  to  the  old  inlet  that  existed  prior  to 
1910.  However,  by  November  1930  the  inner  bar  was  above  mlw  and  the  con¬ 
trolling  depth  over  the  outer  bar  was  2  ft.  By  1935>  the  inlet  was  vir¬ 
tually  closed  and  the  State  of  Delaware  decided  again  to  appeal  to  the 
Federal  Government  for  assistance  in  providing  a  stabilized  inlet. 

As  expressed  by  the  state  and  local  interests,  the  purposes  of  a 
stable  inlet  were  as  follows:  to  provide  for  an  adequate  exchange  of 
water  between  the  bays  and  the  ocean  in  order  to  relieve  stagnation  and 
to  increase  the  salinity  of  the  bays  sufficiently  to  support  a  seafood  in¬ 
dustry,  and  to  provide  e  satisfactory  channel  for  navigation  from  the  ocean 
to  points  within  the  bays.  The  project  designed  by  the  U.  S.  Army  Corps 
of  Engineers  for  these  indicated  purposes  consisted  of  jetties  at  the  in¬ 
let  entrance  extending  to  the  l4-ft  depth  contour  in  the  ocean,  a  channel 
15  ft  deep  and  200  ft  wide  from  the  end  of  the  jetties  to  a  point  7000  ft 
inside  Indian  River  Bay,  and  thence  a  channel  6  ft  deep  and  100  ft  wide  to 
the  natural  6- ft  depth  contour  in  Indian  River  Bay.  This  project  was  ap¬ 
proved  by  Congress,  and  construction  was  completed  in  1939*  A  generally 
satisfactory  inlet  has  existed  since  that  time;  as  of  1966,  a  total  of  ap¬ 
proximately  100,000  cu  yd  of  maintenance  dredging  had  been  performed  and 
the  controlling  depths  in  the  inlet  have  approximated  9  ft,  which  is  ade¬ 
quate  for  navigation.  The  quality  ~f  the  water  within  the  bays  apparently 
has  been  satisfactory  since  opening  of  the  inlet. 

Profile  of  Tide  in  Indian  River  Bay 

One  of  the  assumptions  used  in  the  derivation  of  the  formulas  relat¬ 
ing  to  the  regime  of  tides  in  an  internal  basin  was  that  the  mean  levels 
of  the  water  in  the  basin  and  of  the  sea  outside  are  contained  in  the  same 
horizontal  plane.  It  is  necessary  to  examine  how  nearly  this  assumption 
is  satisfied  for  Indian  River  Bay.  In  table  5,  which  was  prepared  by 
Mr.  Wicker  for  the  guidance  of  the  discussions  herein,  the  longtime  tide 
data  of  the  various  stations  of  the  two  bays  are  given.  The  observations 
of  the  surveys  of  the  two  periods,  1938-1939  and  1948  and  1950,  appear  to 
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be  distinct.  The  former  covers  a  period  of  l4  months.  From  the  1938-1939 
survey  a  superelevation  of  mean  tide  level  (mtl)  of  about  0.15  ft  is  in¬ 
dicated.  The  correspon  dng  value  in  Rehoboth  Bay,  taking  the  average  from 
Dewey  Beach  and  Love  Creek,  is  0.08  ft.  With  reference  to  the  1948  and 
1950  survey  covering  a  period  of  6-1/2  months,  a  much  augmented  super¬ 
elevation  of  mtl  at  Oak  Orchard  and  Dewey  Beach  over  that  observed  in  1938- 
1939  is  indicated. 

The  averaged  1948-1950  values  are  shown  in  fig.  9*  These  were 
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selected  for  detailed  study  in  lieu  of  1938-1939  values  because  current 
velocity  observations  were  made  in  19^8  in  Indian  River  Inlet  and  at  sta¬ 
tions  in  the  bay.  Ocean  values  were  not  observed.  Some  measurements 
would  place  the  range  of  tide  at  the  ocean  as  3*77  ft.  However,  in  a 
memorandum  dated  12  December  1939 >  Mr.  Wicker  rejected  this  value  on  the 
ground  that  the  observations  were  very  poor  and  very  few  in  number.  On 
the  premise  that  the  tide  at  Indian  River  Inlet  is  most  likely  similar  to 
the  tide  at  Lewes,  the  value  4.1  ft  was  proposed  for  the  range.  The  value 
of  tide  at  Lewes  is  well  established.  This  value  is  shown  in  fig.  9j 
H  =  2.05  »  and  will  be  used  in  all  the  subsequent  computations. 

Fig.  9  also  shows  that  the  range  of  tide  along  the  inlet  decreases 
almost  uniformly.  In  the  bay  area,  tidal  displacements  are  nearly  con¬ 
stant.  The  mtl  at  the  Inlet  Bridge  coincides  with  msl,  whereas  in  the  bay 
interior  a  superelevation  of  0.3  ft  is  established.  This  is  a  superele¬ 
vation  much  higher  than  that  observed  during  the  surveys  of  1938-1939* 

The  tributary  discharges,  the  shallowness  of  the  inlet,  and  the  wind  ef¬ 
fects  must  have  a  bearing  on  the  superelevations.  It  is  desirable,  there¬ 
fore,  that  each  of  these  factors  be  separately  examined  in  order  to  es¬ 
tablish  a  reasonable  limiting  value  for  the  superelevation. 


Tribute 


Discharges  into  the  Bays 


In  1930  the  Indian  River  Inlet  Commission  (IRIC)  conducted  a  field 
investigation  to  ascertain  the  inflows  from  the  sea  and  the  outflows  dur¬ 
ing  a  tidal  cycle,  and  also  the  changes  in  the  water  level  within  the 
bays.  The  level  changes  were  noted  for  a  period  of  about  two  years  at 
Holts  Wharf,  and  these  are  shown  in  fig.  10  together  with  the  monthly  pre¬ 
cipitation  in  the  area.  The  flows  were  measured  by  means  of  current 
meters  at  the  Rehoboth  Bay  end  of  Lewes  and  Rehoboth  Canal  and  at  Indian 
River  Inlet.  There  was  no  need  for  flow  observation  in  Assawoman  Canal  as 
this  waterway  was  closed  during  this  study.  Observations  of  currents  were 
made  continuously  each  day  except  Sunday  for  a  period  of  three  months. 
Average  values  from  these  observations  for  consecutive  periods  are  shown 
in  table  6  together  with  the  average  durations  of  inflow  and  outflow. 
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ELEVATION,  FT  RAINFALL,  IN 


PRECIPITATION 


FALL  OF  MEAN  LEVEL  OF  BAY  WATERS 


JUN  JUL  AUG  SEP  OCT  NOV  DEC  JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG 
1929  1930 

Fig.  10.  Bay  water  elevation  and  precipitation 
(Survey  of  Indian  River  Inlet  Commission) 

For  inlet  and  canal  alike,  the  outflow  is  greater  than  the  inflow. 
This  is  attributed  to  the  fact  that  at  the  time  the  average  level  of  the 
water  of  the  bays  was  above  that  of  the  sea.  The  value  of  the  superele¬ 
vation  can  be  determined  analytically,  since  the  disparity  in  the  inflow 
and  the  outflow  durations  is  known.  Let  A  be  the  superelevation,  H  the 
semirange  of  the  tide  in  the  sea,  and  Hg  the  sea  tidal  elevation 
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Fig.  11.  Notation  diagram.  Ocean  tide  and  bay  superelevation 

measured  from  msl,  approximately  H  =  2.05  Ft  (see  fig.  ll).  The  tide  of 
sea  is 


Hg  =  H  sin  ot  ,  o  =  2«/T 

if  t  is  measured  from  the  instant  when  the  tide  is  crossing  the  msl 

plane  and  is  increasing  in  value.  The  period  of  tide  is  T  .  Let  t^  be 

the  first  instant  when  the  waters  of  the  bays  and  of  the  sea  axe  in  the 

same  horizontal  plane.  Let  t.  be  the  duration  of  the  inflow.  Then 
T  1 

t^  +  t^  =  £  t^  is  the  next  instant  when  the  waters  of  the  bays  and  of 
the  sea  are  again  in  the  same  horizontal  plane.  Since  the  internal  tides 
of  the  bays  at  the  time  of  the  flow  measurements  were  negligible,  it  can 
be  written  that 


A  =  H  sin  ot^ 
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and  also  that 


A  =  H  sin  (at1  +  crt^) 


Eliminating  A  between  these  two  equations,  expanding  the  result,  and  col¬ 
lecting  terms,  it  can  be  seen  finally  that 


tan  at^ 


sin  at. 
1 


1  -  cos  at. 

1 


Since  at^  is  a  small  angle 


tan  at^  =  at ^  =  A/&  (68) 

and 

sin  at. 

4*  ■  — 5T5C 

1 

From  table  6  it  is  seen  that  the  average  inflow  duration  is  4.46  hr, 
whereas  the  average  outflow  duration  is  7*l4  hr.  Since  the  sum  of  these 
two  durations,  11.60  hr,  is  less  than  12,42  hr,  the  period  of  the  tides, 
the  observed  inflow  and  outflow  duration  values  need  to  be  increased  1.072 
times.  Thus  the  adjusted  value  of  the  inflow  duration  is  4.79  hr  and  then 
at^  =  139.2  degrees.  Evaluating  the  right  side  of  equation  69  yields 

a/h  =  0.372 


and 


A  =  0.372  X  2.05  =  O.76  ft 


Accordingly,  if  it  is  assumed  that  during  the  study  the  msl  was  1.95  ft 
above  the  1929  datum,  the  mean  level  elevation  in  the  two  bays  should  be 
2.71  ft.  This  is  shown  also  in  fig,  10,  f..  '  it  appears  that  the  value  is 
in  close  agreement  with  the  Holts  Wharf  measurement  of  the  time.  In  the 
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subsequent  analysis  the  notation  Aq  will  stand  for  0 .76  ft 
superelevation . 

Examination  of  the  lower  plot  in  fig.  10  shows  that  the  level  of 
the  bay  waters  fell  about  l/2  ft  from  June  1929  to  April  1930.  It  might 
be  helpful  to  obtain  a  theoretical  relation  to  show  the  dependence  of  the 
rate  of  fall  upon  the  tributary  discharges. 

First  an  expression  is  derived  relating  the  net  total  outflow 
during  a  tide,  for  a  superelevation  A  ,  assuming  in  the  analysis  -chat  the 
variations  of  canal  and  inlet  cross  sections  during  ebb  and  flood  are  ig- 

norable,  and,  more  important,  that  the  tidal  changes  inside  are  insignifi¬ 
cant,  =  A  .  According  to  equation  6  the  inflow  during  a  tidal  cycle 

is  given  by,  after  putting  m  =  1  and  =  A  , 


Q.  =  ka 
1 


%  Aat 
H  '  H  St 


and  the  outflow  by 


where 


Qo  =  ka 


fT+t  1 


£-t 
2  tl 


^2  A 
H~  +  H  dt 


k  -  ■)/  -g.g£*L 

k  ~  V  XL  +  r 

Hp  =  H  sin  ot  ,  a  =  2jt/T 
crt^  =  A/H 


and  a  is  the  cross-/sectional  area  of  canal  or  inlet.  Inspecting  fig.  11 
it  will  be  inferred  that  as  a  good  approximation 

Tr  .  .  1  6  =  t  -  t. 


_ 2  A  _  a\  .  2n0 

h  “  a  \  11 J sin  t 


0  <  9  <  ± 


kh 


and 


and  since  the  square  of  A/H  is  negligible  with  respect  to  unity,  the 
expressions  of  inflow  and  outflow  simplify  to 

and 

%  “r;[1+(i  +  i)i]Tka  <75> 

For  A  positive,  the  inflow  is  smaller  than  the  outflow.  That  is, 

\  -  s.  -  \  (76) 

Introducing  from  equations  7^  and  75  >  the  net  outflow  becomes 

Op  =  ua  TA/fc  (77) 

where 


The  product  pa  is  referred  to  as  the  conveyance  of  an  inlet  or 
canal.  It  is  not  a  dimensionless  quantity;  its  dimensions  are  (pa)  = 

(irT  )  ,  i.e.  volume  per  unit  time.  In  the  early  periods  when  the  Indian 
River  was  being  developed,  the  dimensions  and  hydraulic  characteristics  of 
the  inlet  were  somewhat  uncertain;  therefore,  for  the  analysis  of  the 
tidal  regime  of  that  time,  the  conveyances  could  not  be  evaluated  and  an 
alternate  approach  was  necessary.  During  the  period  from  June  1930  to 
September  1930,  the  net  total  outflow  of  the  inlet  and  canal  was  recorded. 
Let  Qj0  be  this  outflow  and  AQ  represent  the  prevailing  superelevation 
of  the  bay  waters  at  the  time.  From  table  6,  the  combined  value  of  the 
net  outflow  is  Qt0  =  12.76  x  10  ft  .  Analysis  places  AQ  at  0.76  ft. 
Let  be  the  total  net  outflow  at  some  other  superelevation  A  .  Then 


1+6 


and 


Qj  =  T  AQ/fr 


=  Ha  T  A/H 

(it  is  assumed  that  the  longtime  ocean  range  of  tide  is  invariable . ) 

Taking  the  ratio,  it  follows  that 

^  ^  A/Aq  (79) 

o 

showing  that  net  outflow  is  proportional  to  superelevation. 

The  longtime  storage  equation  for  the  two  bays  is,  ignoring  inter¬ 
nal  tides, 

.  dA  _  %  S’ 

dt  T  “  T 

where  is  the  river  inflow  for  one  tidal  cycle.  According  to  weather 

bureau  data,  the  normal  annual  average  evaporation,  as  judged  from  the 
evaporation  in  Class  A  pans,  is  about  39  in.,  whereas  the  normal  annual 
precipitation  is  40  in.  in  northern  Delaware.  Since  these  two  values  are 
nearly  equal,  the  matter  of  precipitation  and  evaporation  does  not  enter 
into  the  storage  equation.  Introducing  from  equation  79, 

A  dA  =  Sr  _  A.  (8o) 

A  dt  T  T  A  v  J 

o 

To  place  the  equation  in  a  dimensionless  form,  first  introduce  T^  ,  the 
number  of  seconds  in  a  year,  and  A^  ,  the  initial  elevation  of  the  bay 
water  at  time  t  =  0  .  Adopting  the  new  variables, 


6  = 

ti  *  Vq 


and  introducing  them  in  equation  80,  the  latter  changes  to 


where 


and 


fq  -  ° i  -  ¥ 


Q,  T 

a  =  •  — 

1  AAj_  T 


B  \  !i  A 

P1  “  A^  '  T  Aq 


(82) 


The  solution  for  6  ,  subject  to  the  condition  that  8=1  when  x  =  0  , 
is 


This  equation  can  be  used  to  determine  the  discharges  from  the  rivers  if 
the  fall  of  the  surface  waters  in  the  bays  is  known  and,  conversely,  the 
fall  of  the  water  surface  if  the  river  discharges  are  known. 

According  to  page  3  of  the  IRIC  report,  the  freshwater  flow  into  the 

6  ^  6 

bays  amounts  to  30  x  10  ftJ  per  day.  This  is  15  x  10°  during  12  hr, 
which  amounts  to  348  ft^/sec.  The  drainage  area  of  the  bays,  including 
their  water-surface  areas,  is  254.5  square  miles.  Deducting  the  water- 
surface  areas,  amounting  to  29.5  square  miles,  assuming  that  evaporation  is 
nearly  of  the  same  value  as  the  precipitation,  the  net  drainage  area  equals 
225  square  miles.  Thus,  the  unit -runoff  is  1.55  ft  per  second  per  square 
mile.  The  only  area  currently  gaged  is  a  tributary  of  Indian  River  with  a 
drainage  area  of  5-24  square  miles.  The  average  discharge  through  1963  is 
741  ftJ/sec,  and  the  unit  runoff  is  1.4l  ftJ  per  second  per  square  mile. 

The  total  drainage  area  of  225  square  miles  does  not  differ  materially  as 
to  runoff  characteristics  from  any  part  thereof.  It  appears  that  the  IR10 
value  o*'  30  x  10°  ft"3  per  day  or  15  x  10  for  12  hr  is  a  reasonably 
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accurate  value  and  will  be  used  in  the  subsequent  computations. 

It  would  be  instructive  to  find  how  closely  the  fall  of  the  water 
level  during  the  period  from  June  1929  to  August  1930  could  be  determined 
using  equation  83.  It  will  be  assumed  that  the  discharges  from  the  tribu¬ 
taries  were  in  proportion  to  the  precipitation  of  the  area,  and  that  during 
this  period  the  conveyances  of  Indian  River  Inlet  and  Lewes  and  Rehoboth 
Canal  remained  constant.  In  particular,  on  the  yearly  average  basis,  0  = 
15  x  10  ft-5  per  tidal  cycle  and  Qrj  =  12.76  X  10°  ft-5  per  tidal  cycle 
when  the  superelevation  is  Aq  =0.76  ft.  Examining  the  upper  graph  in 
fig.  10  it  is  seen  that  the  average  monthly  precipitation  from  June  1929 
to  January  1930  was  2.84  in.,  whereas  from  January  1930  to  August  1930  it 
was  1.64  in.  It  is  now  inferred  that  for  the  former  period  Q  = 

/T  n  zT  o  * 

19.0  x  10°  ft  ,  and  for  the  later  period  =  11.0  X  10°  ft  .  Accord¬ 
ingly,  the  computation  of  the  fall  of  water  surface  in  the  bays  should 
be  made  separately  for  each  period. 

For  the  first  period  it  is  noted  that  the  water-surface  elevation  in 
June  1929  was  about  2.95  ft.  This  makes  A  ^  =  1  ft.  The  relevant  data 
for  the  fall  computation  are: 

A,  =  1  ft  ,  Aq 

O  p 

A  =  8.2  X  10  ft  ,  T/^ 

Or  =  19  X  106  ft3  ,  Qj 

o 

Accordingly,  =  13.8  and  =  15 .1,  and  with  these  values  equation  83 
becomes 


=  0.76  ft 
=  730 

=  12.8  x  106  ft3 


6  =  0.91  +  0.09  e"15’lT 


and  as  A^  =  1  ft 


A  =  0.91  +  0.0 9  e"15'lT 

The  fall  curve  corresponding  to  this  is  shown  in  fig.  10.  In  January  1930 
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the  water-surface  elevation  is  2.86  ft  and  this  makes  A^  =  0.91  ft.  The 
relevant  data  for  the  fall  computation  in  the  second  period  are: 


^  =  0.91  ft  ,  Aq  =  0.76  ft 

O  p 

A  =  8.2  X  10  ft  ,  T/T  =  730 

£  ,  1  (.  ; 

^  =  11  X  10  ff3  ,  Qj  o  12.8  X  10  ftJ 

0 

Accordingly,  0^  =  10.7  and  (3^  =  15-1,  and  with  these  values  equation  83 
becomes 

6  =  0.71  +  29  e“15,lT 


and  as  A^  =  O.91  ft 

A  =  0.65  +  0.26  <f15'lT 

The  fall  curve  corresponding  to  this  is  also  shown  in  fig.  10.  The  agree¬ 
ment  between  the  observed  and  computed  falls  is  sufficiently  close,  sug¬ 
gesting  that  the  use  of  equation  83  is  permissible  to  obtain  the  limiting 
mean  elevation  of  bay  waters  for  the  later  periods  with  the  Indian  River 
Inlet  cross  section  greatly  augmented. 

Limiting  Mean  Elevation  of  Bay  Waters 

With  the  later  enlargement  of  Indian  River  Inlet  the  conveyance  was 
increased  and  the  mean  level  of  the  bay  waters  descended  to  a  lower  plane. 
It  will  now  be  examined  as  to  what  the  limit  would  be  for  the  inlet  con¬ 
dition  present  during  the  surveys  of  1938-1939*  First,  however,  it  is 
necessary  to  establish  the  conveyance  of  the  enlarged  channel. 

According  to  the  IRIC  report,  the  inlet  cross  section  during  the 

g 

study  varied  from  440  to  750  ft  .  Similar  changes  also  occurred  in  the 
depth  of  channel.  For  example ,  the  cut  through  the  barrier  reef  had 
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shoaled  from  -8.0  ft  (zero  at  1929  mlw)  to  which  it  was  dredged  to  an  aver¬ 
age  of  about  -2  ft,  indicating  a  general  fall  of  6  ft.  This  would  mean 
that  the  depth  varied  from  4  to  10  ft.  Taking  the  mean  of  the  extremes, 
the  lower  value  repeated,  the  channel  cross  section  would  he  543  ft  and 
the  depth  6  ft.  It  is  seen,  next,  from  table  7  that  in  1938-1939  the  av¬ 
erage  cross  section  of  the  inlet  channel  in  the  part  next  to  the  sea,  east 
of  sta  22+00,  was  6240  ft  ,  whereas  the  depth  was  11.15  ft.  Gince  the 
conveyance  of  the  channel  would  be  proportional  to 

r2/\ 

then  the  conveyance  of  the  inlet  channel  in  1938-1939  was  about  17.6  times 
as  large  as  the  conveyance  of  the  channel  of  1930.  During  the  period  of 
the  IRIC  study,  the  net  outflow  per  tidal  cycle  through  the  inlet  channel 

6  "3 

as  discussed  previously  was  6.38  x  10  ft  ,  corresponding  to  a  superele¬ 
vation  of  bay  waters  of  O.76  ft  from  msl.  The  corresponding  outflow  capa¬ 
bility  through  the  1938-1939  inlet  channel  would  be  17-6  times  as  large, 
i.e. , 


17.6  x  6.38  X  106  =  112.2  x  106  ft3 

Meanwhile,  since  no  significant  changes  had  occurred  in  the  Lewes  and 
Rehoboth  Canal,  the  net  outflow  per  tidal  cycle  through  the  canal  remained 
the  same  and  thus  it  is 


6.38  x  106  ft3 

Then,  adding  these  two,  the  net  outflow  capability  per  tidal  cycle  in 
1938-1939  through  the  two  waterways  was 

Qj  =  118.6  x  106  ft3,  ^  =  0.76 
o 


The  tributary  discharges  also  remained  the  same: 


/ 


Qr  =  15  X  106  ft3 

To  determine  the  superelevation  of  the  hay  waters,  assume  that  at 
t  =  0  ,  A^  =  Aq  =  O.76.  Also: 

A  =  8.2  x  108  ft2 

=  15  X  108  ft3  per  tidal  cycle 

Qj  =  118.6  x  106  ft3  per  tidal  cycle 
o 

t/pl  =  730 

Hence,  =17-5  and  0^ 

6 

and  since  Aj^  =  0.7 6  ft,  also 

A  =  0.097  +  0.664  6“138t 

Accordingly,  the  terminal  superelevation  which  is  reached  in  about  30  days 
has  the  value  0.097  ft.  This  agrees  fairly  well  with  the  longtime  ob¬ 
served  superelevation  of  the  1938-1939  surveys  at  Oak  Orchard  and  Love 
Creek  (see  table  5). 

Effect  of  Inlet  Cross  Section  on  Bay  Water  Levels 

In  inlets  where  the  mean  depth  of  water  is  small,  the  cross  section 
during  a  flood  would  be  greater  than  that  during  an  ebb.  Assuming  that 
there  are  no  tributaries  leading  to  a  basin  and  that  the  water  level  of 
the  bay  has  attained  its  equilibrium  position,  the  inflow  and  outflow 
volumes  for  a  tidal  cycle  would  be  the  same.  Since  the  currents  of  a 
flood  would  be  greater  than  those  of  an  ebb,  the  duration  of  flood  would 
be  smaller  than  the  duration  of  ebb.  This  circumstance  leads  to  a 


=  138.2,  and  equation  89  reduces  to 
=  0.127  +  0.873  e“138T 
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superelevation  A  which  may  be  related  to  the  half  tidal  range  of  the 
sea  H  .  It  should  be  sufficient  to  examine  this  relation  for  the  simple 
case  of  insignificant  internal  tide. 

Because  surface  width  is  very  large  in  comparison  with  depth,  the 
changes  in  the  cross  section  would  come  from  the  changes  in  the  depth,  r  . 
Thus,  for  the  instantaneous  discharge, 


where  Ai  is  the  differential  head  between  the  bay  and  the  ocean  water 
and  a  factor  of  proportionality  involving  inlet  length,  roughness, 
and  surface  width.  The  depth  of  water  at  mid-inlet  can  be  taken  to  give 
the  effective  hydraulic  radius,  and  on  the  assumption  that  the  surface 
water  in  the  inlet  is  nearly  straight,  its  value  accordingly  would  be 

r  =  rQ  +  ^  sin  0  ,  6  =  at  (85) 

if  the  ocean  tide  is  given  by 


Hg  =  H  sin  ot 

Here  tq  is  the  water  depth  at  time  t  =  0  ,  t  being  measured  from  the 
instant  the  tide  is  rising  from  the  msl  plane.  Equation  85  can  be  written 
in  the  form 

r  =  r  1  +  +  ~~  sin  9 

o  2r  2r 
\  00 

and  since  £/2rQ  is  small,  also 

r  *  ro  (X  +  (X  +  2F  sin  6 )  (86) 

Denoting  by  Ai^  the  instantaneous  discharge  during  an  ebb  and  by  A^ 
the  instantaneous  discharge  during  a  flood,  substituting  in  equation  84 
from  equation  86,  and  recalling  the  meaning  of  H  : 
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(87) 


and 


where 


(88) 


"s  -  «l  V“  ro  (x  +  s; ) 

and  this  is  a  constant. 

T 

At  the  instants  t^  ,  ^  -  t^  ,  and  T  +  ,  the  waters  of  the  bays 

and  of  the  ocean  are  at  the  same  level.  Denote  crt^  by  0^  . 


A  =  H  sin  6^ 


and  since 


0^  is  a  small  angle 


\  =  &/B.  (89) 

Denote  total  inflow  during  an  ebb  by  and  total  outflow  by  Qq  .  From 

equations  87  and  b8: 


and 


Q„ 


sin  0  d0 


5b 


In  the  absence  of  the  tributaries  =  Q,q  and  hence 

Jt-e, 


1  +  2^~  s^n  ®  ysin  d  -  ~  d6 
ro  * 


■  2ji+0. 


it-©. 


1  +  ~~  sin  © 
2ro 


Vt7^ 


sin  ©  d©  (90) 


which  defines  implicitly  the  relation  between  A/R  and  H/2rQ  .  A  simple 
relation  can  be  derived  if  it  is  assumed  that  both  A/R  and  H/2rQ  are 
small  quantities,  so  small  that  the  products  or  the  squares  of  these 
ratios  are  negligible  in  comparison  with  unity.  This  will  be  granted. 
Then  quite  adequately  one  can  put 


Introducing  these  in  equation  90  and  carrying  out  the  integrations  and 
neglecting  in  the  result  the  second-order  quantities,  one  has 


A  1  _H_ 
H  ”5  r 

o 


(91) 


Applying  the  result  to  the  case  of  the  Indian  River  Inlet  with  H  = 
2.05  ft  and  rQ  =  12  ft  one  finds  that 


A  =  0.088  ft 


The  deduction  to  be  made  from  the  discussions  of  the  last  two  sec¬ 
tions  is  that  although  there  are  river  discharges  and  the  inlet  waterway 
is  moderately  shallow,  nevertheless,  the  level  of  the  waters  of  the  two 
bays  is  expected  to  be  only  slightly  elevated.  According  to  the  estimates 
made,  the  superelevation  from  msl  could  hardly  be  greater  than  0.20  ft, 
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a  value  considerably  less  than  that  shown  in  fig.  9,  but  close  to  that 
of  1938-1939?  and  could  very  well  be  ignored  in  the  computations  to  be 
entered  subsequently. 


Wind  Tides  in  Indian  River  Bay 


Examining  again  in  fig.  10  the  course  of  the  elevation  of  the  bay 
waters,  it  is  seen  that  there  was  a  fall  from  September  1929  to  April  1930 
and  then  a  continued  rise  up  to  about  2.80  ft.  This  rise  may  be  explained 
in  part  by  the  shoaling  of  the  inlet.  The  consequent  decreases  in  the  out¬ 
flow  amounts  could  account  for  the  impounding  of  the  waters  of  the  bays. 

Besides  this,  setup  produced  by  the  winds  would  be  another  factor 
to  bear  in  mind.  The  frequently  occurring  strong  winds  can  have  a 
measurable  bearing  on  the  exchange  of  waters  between  the  two  bays  and 
the  ocean.  A  strong  wind  blowing  westward  would  lower  the  water  surface 
at  the  interior  ends  of  inlet  and  canal  and  at  the  same  time  increase  the 
height  of  the  ocean  waters  at  the  exterior  ends  of  the  two  passages.  The 
additional  head  thus  created  may  augment  the  inflow  to  a  degree  that  a  net 
inflow  is  realized  during  a  complete  tidal  cycle. 

Wind  setup  is  governed  by  the  equation 


dA  _  Ts 

dx  ”  Pg  (Dq  +  A) 


T 

S 


=  X  P  V2 
a 


J 


(92) 


Where  A  is  the  setup,  i.e.  the  displacement  of  the  surface  of  water 

from  the  undisturbed  level;  Dq  is  the  depth  of  the  undisturbed  water; 

t  is  the  wind  stress;  p  and  p  are  the  densities  of  water  and  air, 
s  a 

respectively;  V  is  the  wind  velocity;  and  X  is  the  Taylor  stress  coef¬ 
ficient.  The  analysis  of  setup  in  Lake  Erie  due  to  severe  storms  suggests 
that  X  s'  0.0025  (reference  6).  Let  ^  denote  the  lowering  of  waters, 
measured  from  the  undisturbed  level,  at  the  inlet  interior  end.  If  L  is 
the  length  of  the  bays,  the  solution  of  equation  92  subject  to  the  condition 
that  A  is  small  in  comparison  with  Dq  ,  is 
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(93) 


\  *  5  X  p.  vV«  d0p 


Taking  the  average  of  the  lengths  of  the  two  bays  os  the  operative  basin 
lengths ,  L  ■  4.5  miles  (approximately).  Putting  p/pft  -  840  ,  one  finds 
from  equation  93  that  ■  0.045,  0.18,  0.40,  0.72,  and  1.10  ft  corre¬ 
sponding  to  wind  velocities  V  ■  10,  20,  30,  40,  and  50  mpb,  respectively. 

Next,  one  should  consider  the  rise  of  water  surface  at  the  ocean  end 

of  the  inlet  or  canal.  Assume  that  the  Continental  Shelf  is  of  length  L 

and  that  the  depth  of  water  at  the  point  where  the  shelf  terminates  is  of 

value  Dq  .  Assume  that  the  shelf  is  a  uniformly  inclining  surface  of 

inclination  a g  so  that  DQ  ■  .  Since  the  depths  are  great  beyond 

the  shelf,  it  would  be  appropriate  to  suppose  that  wind  setup  vanishes  at 

the  ocean  end  of  the  shelf.  Measure  x  from  the  shoreline,  and  denote 

the  wind  tide  at  as  A  . 

o  o 

Writing 


and 


K 

n 


TsL 


Pg  D, 


XpaV2 


5  -XA  J 


(94) 


the  wind  tide  equation,  equation  92,  takes  on  the  form 


dT)  K 

i 


(95) 


with  the  boundary  conditions 


iBi0  ,  S  -  o 
u  -0  ,  C  « 1 

For  the  problem  at  hand  an  approximate  solution  of  equation  9?  should 
suffice.  As  the  water-surface  height  should  not  differ  substantially 
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from  a  straight  lino  disposition , 


n  «  %(1  •  0 

When  substituted  in  the  right-hand  member  of  equation  95 >  the  result  of  this 
is 


dn  k 

35  ■  ■  10  +  11  •  V  5 

The  solution  of  the  latter,  subject  to  the  condition  that  n  *  0  at 

t  -  1  .  is 


1  -  -j-rV  108  [T|o  +  <X  ’  V  53 

O 


and  as  q  =  q  when  £  =  0  ,  by  substitution 


(96) 


an  expression  that  relates,  implicitly,  the  rise  of  ocean  water  at  the 
coastline  with  the  velocity  of  the  wind. 

Assume  that  the  Continental  Shelf  seaward  of  the  inlet  is  25  miles 
long  and  that  the  slope  is  a  *  l/lOO  .  According  to  equation  9 6  the  wind 
which  should  raise  the  inlet  ocean  waters  by  l/lO  ft  is  of  the  strength  of 
62  mph. 

An  additional  rise  of  water  level  would  be  due  to  the  drift  current 
"rom  the  wind.  The  wind  velocity  being  V  ,  the  drift  current  (i.e.  the 
velocity  of  water  at  the  surface)  would  amount  to  0.03  V  .  Besides  the 
rise  of  the  water  surface  at  the  coastline,  the  wind  will  cause  a  drag  on 
the  waters  of  the  inlet,  augmenting  the  flow  that  would  result  from  a 
difference  of  water  levels  at  the  two  ends  of  the  inlet.  Since  the  com¬ 
putations  needed  to  determine  the  additional  current  are  too  involved, 
they  will  be  ignored  at  this  time. 
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Manning’s  Roughnooo  of  Lewes  and  Rohoboth  Canal 


Tho  main  difficulty  in  the  analytical  evaluation  of  the  surface 
fluctuation  of  a  bay  connected  with  the  ocean  through  an  inlet  or  a  navi¬ 
gation  canal  is  the  assignment  of  a  proper  value  of  the  Manning's  n  for 
the  roughness.  Not  many  references  to  this  question  are  mentioned  in  the 
literature,  and  it  was  thought  to  be  worthwhile  to  attempt  to  obtain  the 
hydraulic  roughness  of  the  Lewes  and  Rehoboth  Canal  on  the  basis  of  the 
field  investigation  results  of  the  summer  of  1930  relative  to  the  amounts 
of  inflow  and  outflow  transpiring  during  complete  tidal  cycles. 

As  indicated  previously,  the  mean  velocity  through  the  canal  can  be 
expressed  as 


where  H  is  half  of  the  tide  range  in  ocean,  r  is  the  hydraulic  radius 
of  the  canal,  and  X  is  the  coefficient  of  friction.  The  connection 
between  X  and  n  is  given  by  equation  15.  The  quantity  is  the 
difference  between  the  levels  of  the  waters  of  Delaware  Bay  and  Rehoboth 
Bay.  Observations  have  shown  that  during  the  period  of  field  investi¬ 
gations  the  internal  tides  of  the  two  bays,  Rehoboth  and  Indian  River, 
were  insignificant.  On  the  other  hand,  the  bay  waters  showed  a  superele¬ 
vation  AQ  with  respect  to  msl.  Analysis  places  this  at  the  value  A0  = 
O.76  ft.  Measuring  from  the  instant  that  the  ocean  tide  is  rising  from 
the  mtl,  during  the  inflow  or  the  flood 

=  H  sin  at  -  A  ,  a  =  2«/T 
and  during  the  outflow  or  the  ebb 

mQ  =  ^  -  H  sin  <rt 

Let  t ^  be  the  first  instant  when  the  waters  of  Rehoboth  Bay  and  of 
Delaware  Bay  are  at  the  same  level.  Then,  the  inflow  during  a  tidal 
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cycle  and  the  outflow 


Introducing  V  from  equation  7  and  denoting  the  canal  cross  section 
during  flood  and  ebb  by  a1  and  a2  ,  respectively,  and  the  correspond¬ 
ing  hydraulic  radii  by  r^  and  rg  ,  one  has 

a  /  2gr  H 

«i  ■  oyurhi  h 

where 


and 


where 


ap  /  2grpH 

%  "  TV  uTZ 


I. 


dot 


The  integrals  1^  and  Iq  were  evaluated  numerically  using  the  data: 
ct^  =0.38  radian,  AQ  »  0.7 6  ft,  and  H  *  2.05  ft.  It  was  found  that 
1±  =  1.44  and  Iq  *  3.45; 


then 


(97) 
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and 


0.29 


(98) 


6  q 

With  the  values  of  the  discharges  known,  i.e.  Q.  =  9*73  x  10  ft-5  and 
Qq  =  l6.1  x  10  ftJ  (see  table  6),  equations  97  and  98  can  be  used  to 
evaluate  two  values  of  X  and  then,  through  equation  15,  the  correspond¬ 
ing  n  values . 

To  complete  the  computation,  the  canal  dimensions  need  to  be  estab¬ 
lished.  The  canal  length  was  about  11.7  miles  between  the  southern  end  of 
the  canal  and  the  northern  end  which  connects  with  Delaware  Bay  via  Broad- 
kill  River  mouth.  At  some  places  the  Lewes  and  Rehoboth  Canal  passes 
through  marshland,  a  circumstance  that  prohibits  an  accurate  description  of 
an  effective  cross  section  from  a  hydraulic  point  of  view.  For  the  problem 
at  hand  it  may  be  sufficient  to  assume  a  uniform  trapezoidal  cross  section 
throughout  the  channel  with  a  bottom  width  bQ  and  with  the  side  embank¬ 
ments  l:m,  and  depth  d  .  In  a  trapezoidal  channel  the  area  equals 
(md  +  bQ)d  ;  the  wetted  perimeter,  2dm  +  bQ  .  It  is  gathered  from  the  IRIC 
report  that  the  bottom  width  was  bQ  =  50  ft  ,  whereas  the  water  depth  was 
6  ft  referred  to  mlw.  Recalling  that  Aq  =  O.76  ft  and  H  =  2.05  ft,  the 
average  depth  in  the  channel  was  d.^  =  9.4  ft  during  high  tide  and 
dg  =  7.4  ft  during  low  tide.  Assuming  that  m  =  2  ,  known  as  a  good 
average  value, 

a]_  =  651  ft2 
rx  =  7.4  ft 

and 

a2  =  480  ft2 

Tg  =  6.0  ft 

Introducing  these  values  in  equations  97  and  98  and  the  proper  values  of 
,  and  L  one  finds  finally  that  during  ebb 

X  =  6.6  x  10-3 
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and  during  flood 


Correspondingly,  during  ebb 
and  during  flood 


X  =  7.4  x  10 


n  =  0.020 


n  =  0.022 


Manning's  Roughness  of  Indian  River  Inlet 

The  velocity  measurements  of  27  July  1948  are  sufficiently  complete 
for  the  determination  of  Manning's  n  of  the  inlet  passage.  The  measure¬ 
ments  were  made  at  the  location  of  the  bridge,  sta  20400.  The  cross  sec¬ 
tion  of  the  inlet  at  this  place  had  been  divided  into  six  parts  and  the 
velocities  in  them  determined  by  current  meters  placed  at  0.6- ft  depths. 
Assuming  that  the  distribution  of  velocities  in  a  vertical  line  is  in  ac¬ 
cordance  with  Manning's  law,  velocity  readings  at  such  depths  would  give 
the  value  of  the  mean  velocities.  The  channel  mean  velocities  as  deter¬ 
mined  from  the  individual  meters  are  entered  in  table  8.  Also  shown  in 
this  table  are  the  time  of  measurement,  the  local  tides,  the  tides  at 
Atlantic  City,  and  the  magnitude  of  instantaneous  cross  sections. 

It  will  be  assumed  that  the  range  of  tide  at  the  ocean  end  of  the 
inlet  is  similar  to  that  at  Atlantic  City.  On  the  day  of  the  velocity 
measurements  the  range  at  Atlantic  City  was  3-30  ft  from  low  to  high  and 
next  2.70  ft  from  high  to  low,  and  3*00  ft  from  low  to  high.  The  mean, 
3.00  ft,  will  be  taken  as  the  apparent  ocean  tide  at  the  inlet  during  the 
day  of  the  observations.  Accordingly,  H  =  1.5  ft.  It  is  preferable  to 
consider  the  quantity  v/  ~\j  2gh  to  study  the  variation  of  mean  velocity 
with  time.  In  this  manner  the  validity  of  the  velocity  data  can  be  ex¬ 
tended  to  any  ocean  tide  range.  Making  the  reductions  of  the  velocity 
data  from  table  8  on  this  basis,  and  using  H  =  1.5  ft,  the  results  are 
shown  in  fig.  12. 

Obviously,  the  velocity  observed  at  the  bridge  should  depend  on  the 
difference  of  the  tides,  Hg  ,  of  the  ocean  at  the  inlet  and  at  the 
bridge.  The  latter  was  observed  and  the  former  may  be  deduced  from  the 
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Atlantic  City  tide  of  the  time  of  the  velocity  measurements.  It  will 

be  assumed  that  the  static  water  levels  in  the  ocean  and  in  the  inlet  at 

the  bridge  are  at  approximately  the  same  elevation  and  that  the  range  of 

tide  at  inlet  ocean  is  nearly  the  same  as  the  range  at  Atlantic  City.  Let 

Hxr  and  denote  the  tides  at  the  inlet  bridge  and  at  Atlantic  CitJ; 

respectively,  when  referred  to  Indian  River  Bay  datum.  These  are  entered 

in  the  second  and  third  columns  of  table  9«  The  average  value  of  H  is 

xr 

2.l6  ft  and  the  average  value  of  Hgr  ,  2.42  ft.  These  differ  but  little 
from  the  corresponding  mtl's  of  2.07  and  2.35  ft.  According  to  the 
assumption,  the  ocean  tide  and  inlet  bridge  tide  with  refer¬ 
ence  to  the  msl  of  the  time  are: 


Hx  =  H 


xr 


2.16 


and 


Hg  =  Hgr  -  2.42 


These  are  entered  in  table  9.  In  the  same  table  are  also  shown  the  values 
of  h  and  h_  ,  representing,  respectively,  the  ratios  H  /H  and 
with  H  =  1.5  ft.  The  variation  of  hx  and  with  time  is 
shown  graphically  in  fig.  13 .  Examining  the  velocity  curve  in  fig.  12, 
it  is  seen  that  the  slacks  occur  approximately  at  the  times  of  hours  9> 
l4.5,  and  21.  At  the  times  of  the  slacks  the  tide  curves  of  ocean  and 
bridge  locale  should  intersect.  This  would  be  fulfilled  if  the  ocean 
curve  based  on  Atlantic  City  tide  in  fig.  13  were  moved  to  the  right  a 
distance  equivalent  to  l/2  hr.  The  curve  thus  adjusted  should  now  repre¬ 
sent  the  ocean  tide  at  the  inlet.  On  this  basis,  the  lag  of  tide  between 
Atlantic  City  and  the  inlet  is  l/2  hr. 


The  velocity  formula  of  equation  7  can  now  be  written  as 
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where  L  is  the  distance  between  the  bridge  and  the  ocean  end  of  the 

X 

inlet,  r  is  the  average  depth  in  the  inlet  between  bridge  and  ocean,  and 
Vffi  is  the  maximum  velocity  in  mid-channel  corresponding  to  the  maximum 
value  of  the  difference  -  h^  .  The  last  equation  can  also  be  written 
as 


From  fig.  12  the  maximum  current  velocities  observed  during  flood,  time 
11.7  hr,  and  during  ebb,  time  17.0  hr,  appear  to  be  of  like  absolute 
values;  that  is, 

Vm 

— -==  =  0.24 
V2gH 

The  differences  in  head  for  the  corresponding  times,  see  fig.  13,  are  also 
of  like  absolute  values;  that  is, 

(*2  -  hx  )  =  °‘27 

\  /  max 

At  time  11. 7  hr  the  cross  section  at  the  bridge,  see  table  8,  is  10,740 

2  2 
ft  ,  and  at  time  18.0  hr,  it  is  9»920  ft  .  Ac  is  expected,  the  mean  of 

2 

these  extremes,  10,330  ft  ,  is  the  same  as  the  cross  section  correspond¬ 
ing  to  low-water  slack  time,  or  the  area  under  msl. 

Since  the  average  cross  section  in  the  inlet  between  the  bridge  and 
the  ocean  is  smaller  than  the  cross  section  at  the  bridge,  the  maximum 
velocity  for  the  inlet  between  the  ocean  and  the  bridge,  a  distance  of 
Lx  =  1430  ft,  -would  be  greater  than  the  value  shown  above.  The  cross- 
sectional  areas  and  the  mean  depths  for  various  years  are  shown  in 
table  7,  and  corresponding  graphs  are  given  in  fig.  14.  From  the  1948 

curve  it  is  estimated  that  the  average  channel  cross  section,  under  msl 

2 

and  between  sta  6+00  and  20+00,  is  8390  ft  .  To  arrive  at  this  value,  the 
stretch  between  the  end  stations  was  divided  into  14  intervals,  the 
areas  at  the  interval  midpoints  were  read  from  the  plot,  and  the  mean 
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was  established.  In  the  same  manner  the  channel  average  depth  was  esti¬ 
mated  to  be  16.8  ft.  The  ratio  of  the  cross-sectional  area  at  the  bridge 
to  the  channel  mean  cross  section  is  1.23.  Thus,  the  effective  maximum 
current  velocity  is 

V 

- —  =  0.24  x  1.23  =  0.296 

Taking  this  and  remembering  that 


max 


0.27 


r  =  16.8  ft 


L  =  1430  ft 
x 


equation  99  yields 


>.  -  246  x  10 


-4 


and  from  equation  15 


n  =  0.046 

As  this  is  a  roughness  value  considerably  higher  than  the  value 
ordinarily  ascribed  to  water  courses,  it  should  be  determined  whether  the 
ocean  tide  range  was  overestimated  when  basing  it  on  the  Atlantic  City 
tide.  The  tide  in  the  ocean  north  of  North  Jetty  at  the  inlet  mouth  3 

observed.  The  observations  as  referred  to  mlw  ocear  datum,  Hgr  ,  are 
entered  in  the  second  column  of  table  10.  Subtracting  2.l6  from  these 
gives  the  ocean  tide  as  referred  to  the  mean  level  of  the  tidal  oscil¬ 
lations  at  the  bridge.  The  value  2.l6  ft  was  established  previously,  see 
table  9-  The  results,  Yl^  ,  are  entered  in  the  third  column  of  table  10. 
The  entries  of  the  last  column  are  h^  ,  the  ratio  H2/fe  ,  H  being  1.5 
ft.  The  data  from  the  last  column  are  plotted  in  fig.  15.  The  curve  of 
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-  OCEAN  TIDE  AT  INLET  BASED  ON 
ATLANTIC  CITY  TIDE 


OBSERVED  TIDE  IN  OCEAN 
NORTH  OF  NORTH  JETTY' 


10  11 
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1,  HR 


Fig.  15.  Comparison  of  inlet  ocean  tides,  observed  and  estimated 


the  ocean  tide  which  was  originally  deduced  from  the  Atlantic  City  data 
is  reproduced  therein  from  fig.  13.  The  two  curves,  one  relating  to 
Atlantic  City  and  the  other  to  the  observations  made  at  the  inlet  mouth, 
are  closely  aligned,  showing  that  determination  from  the  Atlantic  City 
data  is  reliable. 

In  the  above  comparison  the  closer  agreement  is  in  regard  to  the 
range  of  tide  and  not  in  the  temporal  variation  of  tide  for  the  two  deter¬ 
minations.  The  form  of  the  tide  curve  from  the  North  Jetty  observations 
seems  to  be  an  irregular  one  and  therefore  hardly  useful  for  the  present 
computation.  Apparently,  the  determination  of  n  made  in  the  above 
should  be  acceptable.  Since  its  value  is  large  it  would  seem  that  other 
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causes  are  present  to  augment  the  resistance  of  the  inlet  channel.  The 
interferences  of  the  bridge  piers  should  be  considered  if  the  tide  gage 
for  the  bridge  is  placed  west  of  the  bridge.  Fig.  16  shows  that  a  se¬ 
vere  movement  of  sand  is  connected  vrith  the  channel  depth  changes,  and 
these  processes  imply  considerable  energy  dissipation,  perhaps  to  such  a 
degree  as  to  augment  the  resistance  in  the  channel. 

Inertia  Effect  in  Indian  River  Inlet 


In  the  treatment  of  the  problem  of  the  preceding  section  it  was 
assumed  tacitly  that  the  inertia  effect  in  the  flow  through  the  inlet 
channel  is  ignorable.  In  the  event  that  inertia  effect  is  not  ignorable, 
then  at  the  times  of  slacks  there  would  be  a  difference  in  head  between 
the  waters  of  the  inlet  at  the  bridge  and  of  the  ocean.  It  is  now  nec¬ 
essary  to  determine  what  this  diffc.ence  might  be. 

Express  ig  more  precisely,  the  flow  formula  through  a  channel  of 
uniform  cross  section  would  be 


dV 

dt 


dh  V|V| 
"  X  r 


where  the  left  term  is  the  acceleration  of  waters  and  h  is  the  surface 
displacement  measured  in  terms  of  H  ,  the  semirange  of  tide  in  the 
ocean.  Integrating  between  x  =  0  and  x  =  L  : 


L  -h^ 

X  V  V  L 

X 

V  * 

grH 

n  value  at 

slack  times  when  V  =  0  . 

The 

a.  is 

Lx  ( dV  1 

I 

gH  \dt  ; 

t/51* 

'  0 

•  1  (dV) 

(100) 

VsH 

VI#  W  0 

For  the  mean  velocity  of  flow  in  the  channel  at  the  bridge,  fig.  12 
shows  that 
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V 


=  0.2k  sin  at  ,  o  =  2n/T 


y2gH 

The  corresponding  expression  for  the  channel  between  the  bridge  and  the 
ocean  would  oe  1.23  times  as  large,  as  explained  in  the  previous  section, 
and  accordingly, 


U  v  * 

**  V^gH 

Substituting  in  equation  100, 


=  Y  x  0<24  *  1‘23  =  °'59  | 


and  hence, 


^h*)o  =  °-83^ 


L 

Ho  -  H  )  =  0.83  =  —7=  H 

x  'o  T  Yii? 


With  L  =  1430  ft,  H  =  1.5,  and  T  =  4.4?  X  10  sec,  it  follows  that 


(**2  "  Hx)  *  0.018  ft 


0 

which  is  the  maximum  value  for  the  head  fran  inertir  „aring  a  tidal 
period.  It  is  indeed  a  small  quantity  and  can  be  ignored. 


Indian  River  Bay  Tides 


The  formulas  shown  in  the  earlier  sections  were  derived  under  a  set 
of  assumptions,  and  thus  it  is  necessary  to  Know  how  closely  these  assump¬ 
tions  are  verified  in  the  overall  hydraulic  conditions  of  the  bay.  One  of 
the  requirements  of  the  assumptions  was  that  there  are  no  discharges  into 
the  basin  from  tributaries.  There  are  tributaries  emptying  into  the  two 
bays ,  and  by  an  analytical  pi occss  it  was  found  that  the  fall  of  water 

surface  in  1929  and  1930  could  be  eirplained  if  the  volume  of  water  brought 

6  3 

during  a  complete  tidal  period  is  about  15  X  10  ft  .  Since  !-•  e  velocity 
measurements  at  the  inlet  bridge,  surveyed  on  27  July  1948,  indicate  an 
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ebb  outflow  of  nearly  400  x  10°  ftJ  per  tido,  tho  discharges  of  rivers 
into  Indian  River  Buy  can  be  ignored.  A  second  requirement  was  that  the 
mean  .level  of  the  water  in  the  bay  is  in  the  same  plane  as  the  msl.  Hie 
1936-1939  observations  indicate  that  the  longtime  mean  superelevation 
was  0.15  ft  in  Indian  River  Bay  and  0.08  ft  in  Rehoboth  Bay  (see  table  5). 
Analytical  considerations  indicate  that  the  actual  superelevations  may  be 
lower ,  and  it  can  be  assumed  that  there  is  no  significant  superelevation. 

A  final  requirement  is  that  the  connection  of  the  basin  is  with  the  ocean 
through  a  single  inlet  or  with  a  multiple  inlet  system.  This  is  not  true 
for  the  Indian  River  Bay  as  there  is  an  exchange  of  water  between  it  and 
Rehoboth  Bay  through  "The  Ditches."  However,  since  the  volume  of  the 
water  of  the  exchange  is  small  in  comparison  with  the  volume  of  water 
traversing  the  inlet  channel,  one  may  determine  the  Indian  River  Bay  tides 
assuming  first  that  "The  Ditches"  are  closed.  Finally,  these  first  eval¬ 
uations  can  be  corrected  based  on  the  hydraulics  of  "The  Ditches." 

The  determination  of  the  tides  of  the  bay  would  be  in  reference  to 
the  three  quantities:  (l)  the  height  of  mean  high  tide  in  the  bay,  ; 
(2)  the  lag  of  tide,  os  •  and  (3)  the  maximum  current  velocity,  V  ,  in 
the  inlet  passage.  The  formulas  to  evaluate  these  are: 

=  H  sin  t  (51  ) 

a  =  «/2  ”  T  (5:  ) 

\  - 2,1 0 1  f  t  (6:  > 

In  table  4  the  dependence  of  sin  t  and  C  on  the  repletion  coefficient 
K  is  shown.  If,  however,  the  ratio  of  the  inlet  cross  section  a  to  the 
surface  area  of  the  basin  A  is  a  small  quantity,  the  use  of  the  table 
can  be  dispensed  with  since  in  this  case  it  suffices  to  write,  for  K 
less  than  0.3  , 


Thoso  relations  can  be  verified  by  referring  to  table  4  and  fig.  6.  As 
rogards  K  ,  the  coefficient  of  repletion,  it  is  simply,  taking  m  -  1  , 


T  V  2«H  /  r  a 
2«H  *V  XL  +  r  A 


(12) 


the  relation  between  X  and  n  being 


VT 


n  ViT 

1.486  r1/6 


(15) 


As  can  be  seen  from  the  entries  in  table  7,  the  inlet  cross  sections 

are  far  from  uniform.  It  is  necessary  to  normalize  the  cross  sections  to 

arrive  &t  a  prismatic  channel  of  uniform  cross  section,  its  length  the 

same  as  that  of  the  actual  channel,  the  hydraulic  radius  r  equal  to  the 

averaged  hydraulic  radius  of  the  actual  channel,  and  its  n  value  the 

same  as  that  of  the  actual  channel.  Thus,  the  normalization  establishes  a 

new  value  for  the  cross  section.  Using  the  results  of  the  cross-section 

survey  of  1948  it  is  found  that  r  =  12.24  ft  and  a  *»  10,350  ft2  at 

m  m 

msl.  These  values  are  obtained  from  the  1948  curves  in  fig.  14  after  di¬ 
viding  the  stretch  between  sta  6+00  and  60+00  into  20  equal  segments  and 
taking  the  mean  of  the  areas  or  the  depths  for  the  midpoints  of  the  seg¬ 
ments.  Application  of  equation  63  for  normalization  yields  agl  = 

0.944  a_  and  a„0  *=  0.924  a_  for  outflow,  and  thus  the  cross  section 

211  Sc  zu  2 

of  the  normalized  channel  is  a  =  9660  ft  .  Recalling  that  the  n  for 
the  inlet  is  0.046  and  r  *  12.24  ft,  equation  15  yields 


X  -  243  x  10"^ 


The  data  for  the  determination  of  K  from  equation  12  are: 

A  -  420,000,000  ft2 
a  -  9660  ft2 
L  -  5500  ft 
r  -  12.24  ft 
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and 


T  ■  12.42  hr  or  4.47  x  101*  oec 
K  -  2.05  ft 
g  ■  32.2  ft/sec2 
X  -  243  X 

On  the  basis  of  these,  equation  12  gives 

K  -  0.254 

And  using  equation  101,  it  is  found  that 

sin  t  «  0.289 


and 


C  »  0.812 

Theoretically  then,  using  equations  54,  55>  and  62,  the  range  of  tide  in 
the  bay,  2H ^  ,  is  1.19  ft,  the  lag  of  tide  a  =  1.276  radians  or  2.52 
hr,  and  the  maximum  current  velocity  in  the  inlet  is  Vm  =  2.94  ft/sec. 

It  is  seen  from  table  5  that  for  the  years  1948  and  1950  the  mean 
range  of  tide  at  Oak  Orchard  was  0.93  ft  and  at  Maple,  0.99  ft*  The  mean 
of  these,  O.96  ft,  can  be  taken  to  represent  the  range  of  tide  realized 
at  Indian  River  Bay.  This  mean  is  close  to  the  theoretically  computed 
range.  In  table  11  are  shown  longtime  mean  high-  and  low-water  lunitidal 
intervals  for  different  localities.  It  was  argued  in  one  of  the  preced¬ 
ing  sections  that  the  ocean  tide  at  the  inlet  is  about  l/2  hr  later  than 
the  tide  at  Atlantic  City.  Adopting  this  and  since  lag  of  tide  at  Oak 
Orchard  with  reference  to  the  tide  at  Atlantic  City  is  about  3*33  hr,  then 
the  lag  of  the  internal  tide  from  the  ocean  tide  at  the  inlet  should  be 
2.83  hr,  which  is  a  value  not  far  from  the  computed  value.  It  was  shown 
previously  that  the  velocity  survey  of  27  July  1948  gives  for  the  maximum 
current  in  the  bridge  area 
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V  ■  0.24  tJH gli 
m 


With  H  ■  2.05  ft  ,  the  maximum  current  to  bo  had  in  a  cross  section  at 

the  bridge  would  be  V  ■  2.7 6  ft/sec  .  Hinco  the  ratio  of  the  cross 

section  at  the  bridge  to  the  uniform  cross  section  of  the  normalized 

channel  of  the  computations  is  10,330  to  9,660 ,  then  in  the  normalized 

channel  V  "2.95  ft/soc  ,  a  value  which  is  also  close  to  the  computed 
m  ' 

value . 

The  the^  tical  values  of  the  quantities  relating  to  the  internal 
tides  were  corr^id  out  in  the  above,  assuming  that  the  water  .change 
through  "The  Ditches"  is  negligible.  As  there  is  some  exchange  of  water, 
the  next  task  should  be  to  determine  the  effect  of  the  exchange  on  the 
magnitude  of  the  internal  tides. 

The  Hydraulics  of  "The  Ditches" 

Table  11  shows  the  following  data  for  1948: 


Location 

HW  Time,  hr 

Range  of  Tide,  ft 

Inlet  Ocean 

7.54 

4.10 

Oak  Orchard 

10.24 

0.93 

Dewey  Beach 

12.18 

0.47 

It  will  be  assumed  that  the  quantities  observed  at  Oak  Orchard  and  Dewey 
Beach  are  representative  of  Indian  River  Bay  and  Rehoboth  Bay,  respec¬ 
tively.  On  the  basis  of  this  data  the  courses  of  the  tides  in  the  three 
areas  (i.e.  the  ocean,  Indian  River  Bay,  and  Rehoboth  Bay)  would  be  as 
shown  in  fig.  17 •  During  the  entire  time  that  the  level  of  the  water  of 
the  ocean  continues  to  be  higher  than  the  level  of  the  water  of  Indian 
River  Bay,  the  respective  positions  of  water  levels  in  Indian  River  Bay 
and  in  Rehoboth  Bay  undergo  a  reversal.  Thus,  when  an  ebb  of  the  inlet  is 
filling  Indian  River  Bay,  there  first  is  an  inflow  into  Indian  River  Bay 
from  Rehoboth  Bay  and  next,  an  outflow.  The  balance  of  the  discharges  of 
these  two  flows  has  a  bearing  on  the  magnitude  of  the  tides  in  Indian  River 
Bay. 

Let  and  HQ  denote  the  tides  in  Indian  River  Bay  and  in 
Rehoboth  Bay,  respectively,  and  and  the  maximum  values  of 
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Reversal  of  flow  through  "The  Ditches"  during  inlet  flood 


these  tides.  Let  time  t  bo  measured  from  the  instant  when  the  waters 
ol'  Indian  River  Boy  and  the  ocean  are  at  the  some  level  and  the  tide  in 
the  ocean  is  increasing.  Let  tQ  denote  the  log  ol*  tide  in  Rehoboth  Bay 
with  respect  to  the  tide  in  Indian  River  Bay.  The  internal  tides  of  the 
two  bays  ore  very  nearly  sinusoidal.  Actually,  there  are  higher  harmonics, 
but  these  will  be  ignored.  Accordingly,  the  expressions  of  the  internal 
tides  are: 


Hlra  008  0  >  ® 

H  cos  (0  -  t  )  ,  t  »  2*t  /T 
om  '  o  o  o 


(102) 


These  are  the  inferences  from  the  set.  Ignoring  the  slight  dependence  of 
the  lag  of  tide  in  Rehoboth  Bay  with  respect  to  the  tide  in  Indian  River 
Bay  upon  the  magnitudes  of  the  internal  tides,  the  value  of  the  lag  can  be 
placed  at  t  «  2  hr  .  The  ratio  of  the  tides  in  the  two  bays  is  constant, 

“nAm  ’  2  • 

Let  tL  be  the  time  when  the  waters  of  the  two  bays  are  at  the  same 
level.  At  t  =  t^  the  flow  through  "The  Ditches"  reverses  its  direction. 
Place  0-l  =  2nt.j/r  .  Relative  time  0^  is  determined  from  the  condition 


V  li. 


which  in  view  of  equation  102  becomes 

H«m 

COS  0.  =  - —  COS  (0.  -  T  ) 

1  Hlm  10 

Lag  t  was  shown  to  be  2  hr,  and  thus  t  «  60°  .  It  was  also  shown 
o  o 

that  “onA,  ■  l/2  ;  with  this  value  the  last  relation  reduces  to 

cos  0^  ■  cos  0^  +  sin  0^ 
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yielding 


e1  =  60° 

It  is  not  necessary  to  consider  the  flows  through  the  two  "Ditches" 
separately,  and  for  the  purpose  at  hand  one  may  as  well  take  that  the  two 
"Ditches"  can  be  replaced  by  a  single  one  of  cross  section  a  which  is 
the  sum  of  the  cross-sectional  areas  of  the  individual  "Ditc.  ss."  The  ve¬ 
locity  of  inflow  V.  into  Indian  River  Bay  is 


and  the  velocity  of  outflow  is 


Since  a  is  the  cross  section,  the  volume  of  inflow  is 


and  the  volume  of  outflow  is 


Substituting  from  equation  102  and  integrating,  one  has 

Q.  =  0.63  gp  (103) 


and 


Q 

o 


(1CU) 


Accordingly,  the  volume  of  outflow  into  Rehoboth  Bay  during  the  entire 


time  that  Indian  River  Bay  is  being  filled  from  the  ocean  is 

%o  m  °-95  Ir  <105> 

For  the  determination  of  k  ,  which  is  a  quantity  having  dimen¬ 
sions  (l/t-'1'),  proceed  as  follows.  Place 

(106) 

where  CD  is  a  dimensionless  number,  a  cbefficient.  Comparing  the  ex¬ 
pression  for  the  velocity  given  above  with  equation  7  it  will  be 

inferred  that 

CD  (l07) 

The  two  "Ditches"  were  surveyed  in  1935,  and  the  cross-sectional  propor¬ 
tions  are  shown  in  table  12.  Combining  the  two  "Ditches,"  the  average 
value  of  r  when  the  waters  are  at  msl  is  5.8  ft.  The  passages  are  re¬ 
garded  to  be  very  inefficient  for  the  exchange  of  water  between  the  two 
bays.  In  accordance  with  this  idea  it  is  proper  to  assign  to  n  the 
value  of  0.10.  With  r  =  5.8  ft  and  n  =  0.]0  ,  equation  15  yields 
X  =  0.l6l  .  Appearance  of  the  topography  of  the  island  separating  the 
two  "Ditches,"  see  fig.  8,  suggests  that  the  channel  length  of  "The 
Ditches"  is  about  the  same  as  the  maximum  width  of  the  island,  i.e. 

L  =  600  ft.  With  X  =  0.l6l  ,  r  =  5-8  ft  ,  and  L  =  600  ft  ,  equation 
107  yields 


Cp  =  0.24 


and 


k  =  C.24 
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Inserting  this  in  equation  105,  the  volume  of  outflow  into  Rehoboth  Bay 
during  the  entire  time  that  Indian  River  Bay  is  being  filled  from  the 
ocean  is 


«oo—3§§V^  (108> 

This  applies  independent  of  the  magnitude  of  ,  provided  the  tide  in 
Rehoboth  Bay  lags  2  hr  behind  the  tide  in  Indian  River  Bay  and  the  ampli¬ 
tude  of  the  tide  in  Rehoboth  Bay  is  half  the  amplitude  of  the  tide  in 
Indian  River  Bay. 

The  above  treatment  is  introduced  for  the  main  purpose  of  formulat¬ 
ing  the  correction  needed  in  the  evaluation  of  Indian  River  Bay  tides. 
Although  the  flow  through  "The  Ditches"  reverses  its  direction  during  the 
time  Indian  River  Bay  is  being  filled  by  the  flood  of  the  inlet,  for  the 
correction  it-  suffices  to  assume  that  the  flow  is  in  phase  with  the  flow 
through  the  inlet  and  that  the  total  flow  is  of  the  same  magnitude  an 
actually  realized  by  the  currents  of  "The  Ditches."  This  is  permissible. 
Denote  now  the  assumed  hypothetical  flow  through  "The  Ditches"  as: 


This  gives 


or 


V  =  V’  sin  9 
o  om 


Q. 


oo 


sin  9  d© 


Comparing  this  with  equation  108  yields 

=o'ii5V^» 


(109) 
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Indian  River  Bay  Tides  After  Correction 


In  the  presence  of  "The  Ditches,"  the  storage  equation  of  Indian 
River  Bay  is 


dlL 

A“  =  aV-a  V  (110) 

dt  o  o  '  7 


where  is  the  elevation  of  the  waters  of  Indian  River  Bay  referred  to 
msl,  A  is  the  surface  area  of  the  bay,  V  is  the  mean  velocity  of  the 
inlet  currents,  Vq  is  the  mean  velocity  of  the  currents  of  "The  Ditcnes," 
a  is  the  cross-sectional  area  of  the  inlet,  and  aQ  is  the  cross- 
sectional  area  of  "The  Ditches."  From  above 


which,  after  using  the  transformation  relating  to  equations  11,  12,  and 
13,  will  become 


di^ 

w 


> 


h 


The  value  of  the  expressions  contained  in  the  parentheses  is  close  to 
unity,  so  that  in  forming  the  value  it  suffices  to  take  the  velocities 
separately  as  being  sinusoidal.  Consequently, 


V  V’  sin  9 
o  _  om 

V  ”  V  sin  5 

a 


and  the  last  differential  equation  simplifies  to 


-  h  ’  ^  >hl 


(111) 


where 


y;*  \ 

aVa  j 
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The  latter  represents  the  correction  in  the  coefficient  of  repletion  K 

made  necessary  by  the  presence  of  "The  Ditches." 

The  first  estimate  of  H)m  was  =  0.59  ft  *  Using  equation  109, 

V'  =  0.74  ft/sec  .  Also  in  the  first  estimate  it  was  found  that  V  = 
can  '  m 

2.96  ft/sec  .  It  is  reen  from  table  12  that  the  combined  area  of  the  two 

2 

"Ditches,"  surface  waters  at  msl,  is  a  =  5322  ft  .  It  will  he  recalled 
that  the  inlet  cross  section  after  the  normalization  of  the  channel  is 
a  =  9660  ft2.  Hence, 


a  V 

r  r5  -  o-«f> 

m 


and  the  corrected  repletion  coefficient  is 


K’  =  0.86c 


and 


K’  =  0.86  x  0.254  =  0.218 

Accordingly,  one  now  finds  that  the  range  of  tide  in  Indian  River  Bay  is 
1.02  ft  instead  of  1.19  ft  and  the  lag  of  tide  is  2.62  hr  instead  of  2.52 
hr.  These  are  quite  close  to  the  measured  values  of  1948  or  1950.  As 
regards  the  maximum  current  velocity  in  the  inlet,  this  now  is  2.77  ft/sec 
instead  of  2.95,  and  is  an  estimate  lower  than  the  observed  value. 

Erosion  of  Banks  West  of  Highway  Bridge 

Among  the  numerous  manifestations  involving  the  movement  of  sand  in 
Indian  River  Inlet,  such  as  the  formation  of  bars,  depositions  in  the  in¬ 
let  channel,  and  the  stability  of  dredged  forms,  attention  perhaps  should 
be  given  to  the  erosion  of  banks  west  of  the  highway  bridge.  Subsequent 
to  the  dredging  in  June  1938,  it  had  been  observed  that  the  water  widths 
experienced  rapid  and  continued  increases.  The  changes  that  did  occur 
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between  sta  22+00  and  29+00  x'rom  October  1938  to  April  1941  are  shown 
in  fig.  18.  The  underlying  causes  of  this  widening,  a  singular  and  re¬ 
markable  phenomenon,  were  examined  by  Mr.  Wicker.  The  following  exposi¬ 
tion  is  taken  verbatim  from  Mr.  Wicker's  memorandum  of  25  July  1939,  ad¬ 
dressed  to  the  District  Engineer,  Philadelphia  District. 

1.  In  accordance  with  your  instructions,  I  visited  Indian 
River  Inlet  and  made  sufficient  observations  on  July  17,  18,  and 
10,  1939  to  ascertain  the  cause  of  the  erosion  of  Indian  River 
Inlet  channel  banks  west  of  the  highway  bridge.  I  am  convinced 
that  wave  attack  is  responsible. 

2.  The  wav'-s  apparently  occur  only  during  flood  current 
and  the  weaker  parts  of  the  ebb  current.  During  ay  three  days 
of  observation,  the  maximum  wave  height  above  still  water  level 
was  about  one  foot,  the  direction  of  approach  was  approximately 
normal  to  the  bank  throughout  the  embayment  extending  westward 
about  700  feet  from  the  highway  bridge,  and  the  period  of  the 
waves  was  about  15  seconds,  within  +  2  seconds.  These  waves  do 
not  appear  to  have  a  trough,  but  look  and  act  like  mounds  of 
water  moving  rapidly  toward  the  ban!:.  This  appearance,  together 
with  consideration  of  their  period  and  height,  suggests  that  the 
waves  are  of  the  translatory  type.  Translatory  waves  are  often 
derived  waves.  In  this  case,  it  is  believed  that  they  are  de¬ 
rived  from  the  ocean  waves,  which  probably  impart  impulses  to 
the  incoming  flood  current  or  check  briefly  the  feebler  portions 

.  of  the  ebb  current.  When  the  current  is  strongly  ebb,  it  is  be¬ 
lieved  that  the  energy  of  the  ocean  wave  is  rapidly  sapped  by 
the  adverse  current  and  there  remains  no  energy  to  be  dissipated 
by  translatory  waves. 

3.  Current  velocities  and  directions  in  the  embayment 
were  measured.  The  situation  on  flood  current  is,  briefly,  as 
follows:  The  south  embayment  contains  a  counter-clockwise  eddy 
which  reverses  the  flow  along  the  south  bank  from  a  point  about 
4C0  feet  west  of  the  bridge  to  the  bridge.  The  maximum  observed 
current  velocity  in  the  shoreward  side  of  the  eddy  was  about 
7/10  feet  per  second,  which  occurred  at  a  time  when  the  main 
flood  current  thread  had  about  4  feet  per  second  velocity.  It 

is  believed  that  7/10  feet  per  second  is  not  a  scouring  velocity. 
The  embayment  on  the  north  has  a  clockwise  eddy.  Details  were 
not  observed,  but  it  is  believed  that  they  would  have  been  simi¬ 
lar  to  those  found  on  the  south  side.  The  ebb  current  moves 
past  the  embayments  without  creating  any  well-defined  movement 
of  water  within  them.  The  zone  of  moving  water  is  sharply  de¬ 
fined,  as  if  there  were  bulkheads  separating  the  waters  of  the 
embayments  from  the  main  threads  of  current. 

4.  The  waves  do  their  erosive  work  in  the  following  manner: 
On  the  uprush,  they  strike  the  marsh  mud  or  the  scarp  of  the 
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Fig.  l8.  Changes  in  high-water  widths  west  of  Inlet  Bridge 

(Indian  River  Bay) 


spoil  area,  according  to  the  stage  of  tide,  with  considerable 
force.  The  backwash  is  very  strong,  and  sand  moves  steadily 
with  it.  The  marsh  mud  is  attacked,  when  the  tide  is  suffi¬ 
ciently  low,  by  removing  the  sand  layer  under  it.  When  this 
has  occurred,  a  piece  of  marsh  mud  held  together  by  roots  is 
broken  off  and  eventually  disintegrates.  The  material  eroded 
is  deposited  on  the  bottom  in  the  form  of  ridges.  There  were 
three  such  ridges  apparent  a  short  distance  from  the  low  water 
line,  in  depths  of  water  from  0.5  to  1.5  feet.  They  were  about 
four  feet  apart  and  about  6  or  8  inches  high.  Evidently  the 
normal  backwash  cannot  carry  the  sand  that  it  erodes  from  the 
bank  any  farther.  However,  waves  occurring  during  unusually  . 
low  tides,  or  unusually  large  waves,  could  move  the  material 
farther  channelwards ,  where  it  would  be  picked  up  by  the  inlet 
currents . 

5.  Erosion  due  to  wave  attack  can  be  checked  c* ly  by 
longitudinal  works.  In  the  case  at  hand,  revetments  appear*  to 
be  an  adequate  form  of  protection.  There  is  no  evidence  of 
settlement  of  the  rip-rapping  planed  by  the  State,  but  it  is 
being  flanked  on  its  westward  end.  The  work  should  be  done 
as  soon  as  possible,  for  the  erosion  has  now  progressed  to  the 
spoil  bank,  and  there  is  no  reason  for  supp<  sing  it  will  be  less 
rapid  in  the  near*  future .  Delay  will  result  in  the  eventual 
deposition  of  a  considerable  quantity  of  sand  in  the  channel. 

It  would  appear  that  some  of  the  facts  singled  out  in  Mr.  Wicker's 
observations  could  be  explained  on  the  basis  of  turbulent  expansion  of  the 
main  current  entering  into  an  area  with  embayments  on  two  sides.  Ac¬ 
cording  to  analysis,  in  an  expanding  jet  there  will  be  a  lateral  flow  from 
two  sides  induced  by  turbulence  and  directed  to  the  main  current. '  If  the 
embayments  are  limited  in  area,  the  longitudinal  flow  at  the  edges  of  the 
expanding  jet  together  with  the  lateral  flow  to  the  edges  will  enhance  the 
formation  of  eddies  in  each  embayment.  The  eddies  were  observed. 

In  a  current  moving  between  two  embayments  and  experiencing  lateral 
turbulent  expansion,  the  current  velocities  decrease  when  moving  laterally 
from  the  center  line  of  the  current  toward  embayment  banks.  Also,  very 
likely  the  depth  of  water  decreases  from  the  center  line  away  to  the  banks. 
These  two  conditions  will  cause  the  refraction  of  waves.  If  U  is  the 
current  velocity  and  d  the  depth  at  a  place,  a  point  on  the  wave  crest 
would  be  moving,  in  space,  with  velocity 

w  =  s/9d  +  U 
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Accordingly,  when  a  solitary  wave  with  crest  originally  normal  to  the 
inlet  channel  walls  enters  into  the  area  with  embayments,  the  point  of 
the  crest  at  the  center  moves  away  with  greater  velocity  of  propagation 
and  points  farther  away  from  the  center  with  lesser  velocities.  Conse¬ 
quently,  the  original  wave  breaks  into  two  waves  eventually,  with  align¬ 
ments  parallel  to  the  banks  and  moving  toward  the  bank.  This  also  was 
observed. 
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Table  2 

Tabular  Values  of  the  Coefficient  of  Repletion  for  T  =  12  hr 
and  n  =  0.02,  0.03.  0.04,  and  0.05 


_ (AT,  VH/a)  x  10-4 _ 

X,  =  0.1  X  g  L  =  0.2  X  IQ1*  L  =  0.5  X  10* L  s  1.0  X  10*  L  =  2.0  X  10*  L  =  4.0  X  IQ* 

n  =  0.02 


5 

3.697 

2.859 

1.974 

1.443 

1.11? 

0.741 

10 

4.557 

3.825 

2.867 

2.180 

1.606 

1.160 

15 

4.810 

4.364 

3*437 

2.708 

2.042 

1.495 

20 

5.004 

4.611 

3-830 

3.U0 

2.397 

1.780 

30 

5-293 

5.018 

4.637 

3.722 

2.952 

2.254 

4o 

5.295 

5.120 

4.618 

4.053 

3.209 

2.626 

60 

5.383 

5-327 

4.938 

4.537 

3.909 

3.195 

8o 

5.426 

5-338 

5.101 

4.770 

4.262 

3-596 

100 

5. >*49 

5.383 

5.195 

4.525 

4.500 

3.899 

n  =  0.02 


5 

2.734 

2.064 

1.365 

O.98O 

O.699 

0.496 

10 

3-704 

?  975 

2.098 

1.519 

1.097 

0.783 

15 

4.219 

3-^3 

2.589 

1.939 

1.4i6 

I.019 

20 

4.528 

3-933 

2.984 

2.284 

1.689 

1.223 

30 

4.871 

4.4i4 

3-556 

2.825 

2.29b 

1.577 

40 

5.051 

4.688 

3-941 

3.230 

2.508 

1.873 

60 

5.233 

4.985 

4.422 

3.794 

3-067 

2-358 

80 

5-323 

5-140 

4.696 

4.154 

3.473 

2.743 

100 

5-366 

5.228 

4.870 

4.406 

3.778 

3.055 

1 


n  =  0.04 


5 

2.171 

1.598 

1.038 

0?4o 

0.525 

0.385 

10 

3-101 

2.389 

1.604 

1.158 

0.828 

0.589 

15 

3-732 

3-264 

2.040 

1.495 

1.077 

0.769 

20 

4.044 

3.344 

2-397 

1.780 

1.292 

0.927 

30 

4.504 

3.900 

2.948 

2.245 

1.662 

1.134 

40 

4.761 

4.252 

3-328 

2.626 

1-971 

1.441 

60 

5-037 

4.665 

3.909 

3.191 

2.473 

1.846 

80 

5.175 

'  4.894 

4.262 

3.598 

2.866 

2.180 

100 

5-261 

5-035 

4.500 

3.894 

3.182 

2.460 

5 

1.793 

1.3U 

0.836 

0.594 

0.422 

0.298 

10 

2.675 

1.979 

1.303 

0.935 

0.666 

0.473 

15 

3.202 

2.668 

1.675 

1.213 

0.869 

0.637 

20 

3-603 

2-873 

1.985 

1.452 

1.036 

0.74y 

30 

4.133 

3-447 

1.580 

1.858 

1-353 

0.973 

40 

4.454 

3-837 

2. 8*1 

2.194 

1.615 

1.168 

60 

4.819 

4.336 

3.452 

2.726 

2.055 

1.512 

80 

5.012 

4.622 

3-843 

3.128 

2.412 

1.796 

100 

5-123 

4.8*5 

4.331 

3-392 

2.713 

2.046 

Table  3 


Parameters  of  the  Fluctuation  of  Water  Surface 


in  the  Basin  as  a  Function  of 

_K 

K 

al 

_a3_._ 

b3 

COS  T 

sin  t 

0.1 

0.9936 

-0.0001 

-0.0052 

0.99327 

0.115804 

0.2 

0.9745 

-0.0004 

-0.0106 

0.97334 

0.22934 

0-3 

0.91+35 

-0.0009 

-0.0164 

0.94086 

0.33874 

0.U 

0.9020 

-0.0017 

-0.0220 

0.89735 

0.44137 

0.5 

0.8515 

-0.0028 

-0.0282 

0.84425 

0.53593 

0.6 

0.7942 

-0.0043 

-0.0347 

0.78386 

0.62091 

o.? 

0.7325 

-0.0063 

-0.0418 

0.71856 

0.69549 

0.8 

0.6689 

-0.0089 

-0.0495 

0.65091 

0.75917 

0-9 

0.5997 

-0.0123 

-0.0579 

0.57732 

0.81649 

1.0 

0.5451 

-0.0165 

-0.0664 

0.51783 

0.85551 

1.2 

0.4369 

-0.0281 

-0.0849 

0.39949 

0.91676 

1.4 

0.3489 

-0.0448 

-0.1038 

0.30119 

0.95357 

1.6 

0.2811 

-0.0661 

-0.1201 

0.22449 

0.97446 

1.8 

0.2294 

-0.0910 

-0.1327 

0.16588 

0.98614 

2.0 

0.1893 

-0.1177 

-0.1401 

0.12160 

0.99258 

3-0 

0.8830  x  10*1 

-0.2207 

-0.1187 

0.02953 

0.99956 

4.0 

0.5032  X  10-1 

-0.2606 

-0.0802 

0.01037 

0.99995 

5-0 

0.3232  x  10"1 

-0.2740 

-0.0532 

0.00575 

0.99898 

6.0 

0.2249  x  10- 1 

-0.2794 

-0.0377 

0.00363 

0.99999 

7.0 

0.1653  x  10"1 

-0.2817 

-0.0280 

0.00256 

1.0000 

8.0 

0.1266  X  10'1 

-0.2828 

-0.0215 

0.00192 

1.0000 

9-0 

0.1001  x  io_ 1 

-0.2835 

-0.0170 

0.00150 

1.0000 

10 

0.8105  X  10-2 

-0.2845 

-0.0138 

0.00119 

1.0000 

20 

0.2026  x  10"2 

-0.2845 

-0.0035 

0.00030 

1.0000 

30 

0.9007  X  10-3 

-0.2845 

-0.0015 

0.00013 

1.0000 

4o 

0.5066  x  10-3 

-0.2845 

-0.0009 

0.00008 

1.0000 

50 

0.3242  x  10-3 

-0.2845 

-0.0006 

0.00005 

1.0000 

6o 

0.2252  x  10-3 

-0.2845 

-0.0004 

0.00004 

1.0000 

70 

0.1654  X  10-3 

-0.2845 

-0.0003 

0.00003 

1.0000 

80 

0.1267  X  10-3 

-0.2845 

-0.0002 

0.00002 

1.0000 

90 

0.1001  x  10-3 

-0.2845 

-0.0002 

0.00001 

1.0000 

100 

0.8105  X  10-4 

-0.2845 

-0.0001 

0.00001 

1.0000 

Table  4 

Coefficient  C  in  Tidal  Prism  Formula 
and  Range  of  Tide  in  Basin 


K 

sin  t 

C 

K 

sin  t 

c 

0.1 

0.1158 

0.8106 

4.0 

0.9999 

0.9993 

0.2 

0.2293 

0.8116 

5.0 

0.9999 

0.9994 

0.3 

0.3387 

0.8128 

6.0 

1.0000 

0.9997 

0.4 

0.4414 

0.8153 

7.0 

1.0000 

0.9997 

0.5 

0.5359 

0.8184 

8.0 

1.0000 

0.9998 

0.6 

0.6209 

0.8225 

9.0 

1.0000 

0.9998 

0.7 

0.6955 

O.8288 

10.0 

1.0000 

0.9998 

0.8 

0.7592 

0.8344 

20 

1.0000 

0.9998 

0.9 

0.8165 

0.8427 

30 

1.0000 

0.9999 

1.0 

0.8555 

0.8522 

4o 

1.0000 

0.9999 

1.2 

0.9168 

0.8751 

50 

1.0000 

0.9999 

1.4 

0.9536 

0.9016 

60 

1.0000 

1.0000 

1.6 

0.9745 

0.9267 

70 

1.0000 

1.0000 

1.8 

0.9861 

0.9484 

80 

1.0000 

1.0000 

2.0 

0.9926 

0.9650 

90 

1.0000 

1.0000 

3-0 

0.9996 

0.9950 

100 

1.0000 

1.0000 

Table  5 

Longtime  Mean  Data  of  Tides  in  Indian  River  and  Rehoboth  Ba; 
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Table  8 

Indian  River  Inlet  Velocities 
Survey  of  27  July  1948 


Time,  hr 
e.s.t. 

Atlantic  C 

Tide 

ity  Inlet  Bridge 

Velocity 

ft/sec 

Cross  Section 
ft2 

Discharge 

ft^/sec 

8.95* 

2.60 

2.14 

0.00 

10,360 

0 

9.91 

3.00 

2.41 

+1.29 

10,470 

+13,510 

10.0? 

3.32 

2.66 

+1  . 

10,580 

+20,100 

10.64 

3.53 

2.83 

+2.16 

10,660 

+23,030 

11.20 

3-69 

2.95 

+2.30 

10,710 

+24,630 

11.76 

3.70 

3.00 

+2.36 

10,740 

+25,350 

12.32 

3.56 

2.98 

+2.17 

10,730 

+23,280 

3.32 

2.91 

+1.84 

10,700 

+3.9,690 

13.4p 

3.01 

2.82 

+1.43 

10,660 

+15,240 

14.01 

2.69 

2.68 

+0.87 

10,590 

+9,215  . 

14.57** 

2.30 

2.44 

0.00 

10,490 

0 

15.20 

1.90 

2.12 

-1.17 

10,350 

-12,110 

15-83 

1.50 

1.80 

-1.90 

10,210 

-19,400 

16.45 

1.20 

1.51 

-2.23 

10,08C 

-22,480 

17.08 

1.01 

1.28 

-2.36 

9,975 

-23,540 

17.71 

1.04 

1.15 

-2.30 

9,920 

-22,820 

18.34 

1.21 

1.19 

-2.10 

9,940 

-20,870 

18.97 

1.47 

1.31 

-1.82 

9,990 

-18,180 

19  59 

l.f: 

1.50 

-1.45 

10,070 

-14,600 

20.22 

2.24 

1.75 

-0.80 

10,180 

-8,145 

20.85* 

2.68 

2.04 

0.00 

10,310 

0 

*  Low-\/ater  slack. 
**  High-water  slack. 
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Table  11 

Longtime  Averages;  High-  and  Low-Water  Times  and  Ranges 


Location _ 


June,  Julv 


Time  of 
High  Water 
hr,  e.s.t. 


Time  of 
Low  Water 
hr,  e.s.t. 


Range  of 
Tide 
ft 


and  September,  1948  Survey 


Atlantic  City 

7-24 

1.12 

4.08 

Lewes 

8.46 

2.11 

4.14 

Indian  River  Inlet  Bridge 

7-7  6 

1.47 

2.55 

Oak  Orchard 

10.24 

4.65 

:  0.93 

Dewey  Beach 

12.18 

6.43 

0.47 

January,  February,  March,  and  April,  1950  Survey 


Atlantic  City 

7.24 

1.12 

4.08 

Indian  River  Inlet  Bridge 

7.89 

1.60 

2.60 

Oak  Orchard 

10.53 

4.63 

0.93 

Maple 

10.93 

5-50 

0.99 

Table  12 

Cross-Sectional  Dimensions  of  "The  Ditches" 


5  S 


Water  Level* 

2  Su 

rfaee  Width 

Mean  Depth 

ft 

Area,  ft 

ft 

ft 

Little  Ditch 

2.8 

3150 

420 

7.5 

2.0  (msl) 

2816 

4l6 

6.8 

0.0 

2000 

Big  Ditch 

380 

5.3 

2.8 

3620 

900 

4.0  I 

2.0  (msl) 

2506 

894 

2.8  I 

0.0 

1114 

85G 

*  Level  referred  to  1929  datum. 
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